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Passive red galaxies frequently contain warm ionized gas and have spectra similar to low-ionization 
nuclear emission-line regions (LINERs). Here we investigate the nature of the ionizing sources pow- 
ering this emission, by comparing nuclear spectroscopy from the Palomar survey with larger aperture 
data from the Sloan Digital Sky Survey. We find the line emission in the majority of passive red 
galaxies is spatially extended; the Ha surface brightness profile depends on radius r as r -128 . We 
detect strong line ratio gradients with radius in [N Il]/Ha, [S Il]/Ha, and [O III] / [S II], requiring the 
ionization parameter to increase outwards. Combined with a realistic gas density profile, this out- 
ward increasing ionization parameter convincingly rules out AGN as the dominant ionizing source, 
and strongly favors distributed ionizing sources. Sources that follow the stellar density profile can 
additionally reproduce the observed luminosity-dependence of the line ratio gradient. Post-AGB stars 
provide a natural ionization source candidate, though they have an ionization parameter deficit. Ve- 
locity width differences among different emission lines disfavor shocks as the dominant ionization 
mechanism, and suggest that the interstellar medium in these galaxies contains multiple components. 
We conclude that the line emission in most LINER-like galaxies found in large aperture (> 100 pc) 
spectroscopy is not primarily powered by AGN activity and thus does not trace the AGN bolometric 
luminosity. However, they can be used to trace warm gas in these red galaxies. 

Subject headings: galaxies: active — galaxies: ISM — galaxies: elliptical and lenticular — galaxies: 
emission lines — ISM: kinematics and dynamics — stars: AGB and post-AGB 



1. INTRODUCTION 

Emission lines are important spectral features that can 
help us probe the gaseous component in galaxies. They 
are not unique to star-forming galaxies, but also exist 
in galaxies with only old stellar populations. Numer- 
ous result s (iPhillips et all 119861 : IGoudfrooii et all 119941 : 
I Yan et all 120061 ) have shown that line emission is preva- 
lent in more than 50% of passive red galaxies, and they 
have line ratios similar to t he low ionization nuclear emis- 
sion line regions (LINERs, EeckmadEHn). What pow- 
ers this line emission has been an unsettled question for 
decades. 

LINERs are identified by their particular pattern 
of line strength ratios, with strong low-ionization for- 
bidden lines (e.g. [Nil], [S II], [Oil], [Ol]) rela- 
tive to recombination lines (e.g. Ha, H/3) and high- 
ionization forbidden lines (e.g. [O III]). Unlike galax- 
ies dominated by star-forming HII regions or Seyfcrts, 
whose line ratio patterns clearly identify their ioniz- 
ing sources as young massive stars or active galactic 
nuclei (AGN), respectively, LINERs can be produced 
by a wide array of ionizatio n mechanisms, such as 
photo-ionization by an AGN dFerland fc Netzerl 119831: 
lHalpern fc Steinerlll98l iGroves et al.ll2004bD. p hotoion- 
ization by post- AGB stars dBinette et alJll994[ h fast ra- 
diative shocks (jDopita fc Sutherland! I1995D. photoion- 
ization by the hot X-ray-em itting gas ()Voit fc Donahue! 
ll99CtlDonahue fc Voitlll991f) . or thermal conduction from 
the hot gas ( jSparks et alj |1989). Therefore, their exact 
ionization mechanism has been hotly debated. 

The LINER puzzle is futher complicated by the lim- 
ited spatial resolution available in many samples, par- 



ticularly in SDSS. Originally, "LINER" only referred 
to a class of galaxy nuclei. They we re first identifie d 
in nuclear spectra of nearby galaxies (|Heckmanl 11980) . 
This is t he case in most L INER studies of nearby galax- 
ies fe.g. IHo et all 11997b!) . For our discussion, we re- 
fer to these LINERs as "nuclear LINERs." We should 
keep in mind that ground-based slit spectra, under typ- 
ical seeing, usually cannot resolve better than the cen- 
tral few hundred parsecs for even nearby galaxies (< 40 
Mpc). This is the scale referred to by the word 'nuclear'. 
With na rrow band ima g ing and / or l o ng-sli t spectroscopy 
survey s iPhillips et al.l (fl986h. iKiml (U989T). iBuson et al.l 
(fl993fl. IGoudfrooii et all (|19940Macchetto et al.l(|1996lL 



iZeilinger et al.l (|1996l ). and others found that the line 
emission in early-type galaxies is often extended to 
kpc scale and has LINER-like line ratios. We refer 
to these cases as "extended LINERs." In surveys of 
much more distant galaxies, such as most galaxies in 
SDSS, or surveys at high red shifts, such as DEEP2 
([Davis et al.ll2003lV zC OSMOS (jLillv et al.ll2007f l. BOSS 
(lEisenstein et al ]l 2011f) . and distant cluster surveys (e.g. 
Lubi n et al.l 120 09). the spectra obtained usually covers 
a much larger scale than the nuclei and we are not al- 
ways able to tell how the emission is distributed spatially. 
Still, a large num b er of galaxies show L I NER-like spectra 
(lYan et al.l 120061: iLemaux et al.1 120101: iBongiorno et al.l 
l2010HYan et al.ll2011l ). The term "LINER" is often casu- 
ally adopted in this case to refer to galaxies with LINER- 
like spectra. Here, we refer to these cases as "LINER-like 
galaxies." 

Although the name of LINER includes a morphological 
description — "nuclear" , there has been no quantitative 
definition of what line emission distribution would qual- 
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ify as a nuclear LINER. All LINERs have been defined 
only spectroscopically based on their line ratio patten. 
The distinction between nuclear LINERs and extended 
LINERs is very murky. Practically, it often refers to the 
scale over which the spectrum is taken — < 200 pc for nu- 
clear and > 1 kpc for extended, rather than a characteris- 
tic scale or m orphological des c riptio n of the line emission 
distribution. iMasegosa et al.l (|2011[ ) tried to classify the 
different morphologies of nuclear line emission distribu- 
tion based on narrow band HST images. In the large 
majority of LINERs (84%), they found an unresolved 
nuclear source surrounded by diffuse emission extending 
to a few hundreds of parsecs. However, it is unclear in 
those data whether the two components have different 
ionizing sources and which component dominates in to- 
tal luminosity. 

For a large fraction of nuclear LINERs, which are de- 
fined using spectra from the central few hundred par- 
secs, AGN activity clearly exists in the center, although 
some puzzles still remain. Evidence for AGN activity in- 
cludes the detection of central hard X-ray point sources, 
compact radio cores, broad emission lines in d irect or 
polarized light, and UV variability ( see lHdl2008l and ref- 
erences t herein). On the other hand, there is growing 
evidence (jHo et al.ll20"0lt lEracleous et al.ll2010D that the 
weak AGN in most nuclear LINERs does not emit enough 
photoionizing photons to account for the observed in- 
tensity in optical emission lines. Based on narrowband 
images or slit spatial profiles observed from HST, the 
narrow line region in LINERs appear to be strongly con- 
centrated in the cente r, with typical dim ensions smaller 
than tens of parsecs ()Walsh et al.l [2008) . However, in 
some o bjects the profile c an extend to a few hundred 
parsecs ( Shi elds et al.l[200l . Thus, even for nuclear LIN- 
ERs, whether the AGN is responsible for all of the narrow 
line emission within the central few hundred parsecs is 
unclear. 

For the extended LINERs, it is totally unclear how 
they are related with nuclear LINERs and what mecha- 
nism powers their line emission. For LINER-like galax- 
ies, we know even less. Are they very powerful nu- 
clear LINERs or are they dominated by extended LINER 
emission? The nuclear-LINER fraction among early- 
type galaxies is fairly similar to the LI NER-like galaxy 
fraction among red galaxies in SDSS (|Ho et al.l Il997ri 
I Yan et all 120061 ) . Is this similarity fortuitous, or is there 
a physical connection? We try to answer these questions 
in this paper. 

Solving these puzzles is important. Many people have 
used the line strength in LINER-like g alaxies from SPSS 
as a n indicator of AGN p o wer ( I Kauffmann et al 
20031: iConstantin fc Vogelejl 12009: iKewlev et al' 
2006 : iSchawinski et al.l l2007t IKauffmann fc Heckman 
2009f ). while others ar gued they are not g enuine 
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but more likely powered by hot evolved stars. Therefore, 
settling this puzzle is not only important for AGN 
demographics, but also for understanding the ISM in 
early-type galaxies. If the line strength is not an AGN 
indicator, it might instead reflect the amount of cool gas 
supply, or the cooling rate of the hot gas. To make the 
correct physical connections, we have to first find out 
what powers the line emission. Additionally, if we can 



pin down the ionization mechanism, we will be able to 
measure the gas-phase metallicity in these red galaxies, 
which has so far been impossible. 

For the extended LINER emission, three types of ev- 
idence have been presented to argue that they are not 
powered by AGN. (a) Po st-AGBs could produce enough 
photo-ionizing pho t ons (Idi Serego Alighieri et al.l 119901 : 
IBinette et al.j|1994t iStasihska et al.ll2008D . (b) The line 
emission regions are spatially extended wit h a surface 
brigh tness profile that is shallower than r~ 2 (jSarzi et al.l 
I2010T) . (c) The line luminosity correlates with ste llar lu- 
minosity (|Sarzi et al.ll20Tol : ICapetti fc Baldill20ll . This 
evidence led to the reasonable suspicion that a popula- 
tion of hot evolved stars might in fact dominate over the 
AGN photoionization in early-type galaxies. However, 
all of these arguments depends on inherent assumptions 
about observationally unknown factors. 

For the first argumen t, Idi Serego Alighieri et all 
( 1990). IBinette et al.l (|1994f ). and IStasihska et al.l (|2008f ) 
all assumed that nearly all the photons produced by post- 
AGB stars are absorbed by the gas and that the gas is 
distributed near the ionizing stars to produce the correct 
ionization parameter. Both the gas covering factor and 
the relative distribution of gas to stars are unknown. 

For the second argument, ISarzi et ail (|2006l 12010ft 
showed that the line emission in nearly all line-emitting 
early-type galaxies is spatially extended, and has a sur- 



face bright n ess p rofile shallower than 
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ISarzi et all JloTfl) is very careful in not drawing con- 
clusions based on this fact alone, such e xtended emis- 
sion has been quoted by many others (e.g. rKaviraill2010l : 
iMasters et~all 120101: ISchawinski et al.lboiOl ) as evidence 
for stellar photoionization. However, a central point 
source can also produce extended line emission regions. 
The line emission surface brightness profile not only de- 
pends on the ionizing flux profile, but also depends on 
how the gas filling factor, spatial distribution of gas 
clouds, and the gas density vary with radius. Both the 
filling factor and the cloud spatial distribution are so 
poorly known that the surface brightness profile provides 
no practical constraint on the flux profile. 

The third argument is based on the observed corre- 
lation in surface b rightness between line emission and 
stellar continuum. Sarzi et al. (2010) showed with IFU 
data that the H/3 equivalent width (EW) is fairly con- 
stant th roughout each ga l axy w hen excluding the nuclear 
region; ICapetti fc Baldil (|2011l ) argued a similar point 
based on integrated line emission from SDSS. This might 
seem like the strongest support for stellar photoioniza- 
tion. However, as mentioned above, the line flux de- 
pends on many other unconstrained factors besides the 
ionizin g flux, such as the gas filing factor. ISarzi et all 
(2010) performed the calculation for a simple model of 
stellar photoionization and found that it did not produce 
the spatially-invariant EW of H/3 they found through- 
out each galaxy in their sample. To make the model 
consistent with observations, certain fine tuning of the 
gas filling factor, density, and/or mean-free path of the 
ionizing photons is required, which has no independent 
observational support. 

In fact, as long as both the line flux and stellar contin- 
uum are strong functions with radius in these galaxies, 
one would always find tight correlation between the two, 
even in the case of photoionization by an AGN. 
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Thus, we need a simpler test that can distinguish 
different ionization mechanisms that is free of assump- 
tions about unknown parameters. A central point source 
(e.g. AGN) and a system of spatially distributed ionizing 
sources will produce different ionizing flux profiles as a 
function of radius. Thus, for the same gas density profile, 
they would yield different ionization parameter profiles, 
leading to different spatial gradients in line ratios. Thus, 
looking at the line ratio gradients may provide a clue 
to differentiate the two scenarios. The only assumption 
involved here is the gas density profile, which can be de- 
rived from the hot gas density profile assuming pressure 
equilibrium or can be measured directly from line ratios. 

Ideally, the line ratio gradient is best measured from 
integral field spectroscop y (IFS) data. Ho wever, cur- 
rent IFS data (SAURO N, Bacon et al.l l200lt ATLAS3D, 
Cappcllari ct al. 20TT1) have very limited wavelength cov- 
erage and do not probe enough strong emission lines 
to detect the ionization gradient in a large number of 
galaxies. The emission lines covered by SAURON are 
O III] , H/3, and in a few cases [N I] A5197,5200. Because 
O III] /H/3 depends on the ionization parameter and to 
some extent also on the hardness of the ionizing spectra, 
it alone does not provide an unambiguous constraint on 
the ionization gradient. In addition, the weakness of H/3 
makes it more difficult to detect small changes. 

Our solution i s to use the nu clear spectra from the 
Palomar survey ()Ho et al.| [T995). and the fiber spectra 
from SDSS. By identifying the same population of galax- 
ies at different redshifts, for which the fixed angular aper- 
ture corresponds to different physical radii, we can study 
statistically the spatial profile of emission line surface 
brightness and the line ratio gradient. The emission line 
surface brightness profile can also inform us about the 
relationship between nuclear LINERs and extended LIN- 
ERs. 

In addition, we will examine the widths of emission 
lines. Different forbidden lines will have different widths 
if there is line ratio variation within a galaxy and the 
variation is correlated with gas kinematics. On the other 
hand, shock ionization models also should produce width 
differences among multiple lines due to the dependence 
of line ratios on shock velocity. 

Our investigation includes all line-emitting red galax- 
ies except for those containing star formation; we do not 
specifically select for LINER-like galaxies. Based on line 
ratio diagnostics, these line-emitting red galaxies do have 
fairly uniform line ratios with the most typical case be- 
longing to the LINER category, as shown by the line ratio 
diagnostic diagram in Figure [TJ Avoiding the use of line 
ratios in sample selection is essential for our study of the 
line ratio gradient. 

This paper is organized as following. In Section 2, we 
will describe the data, measurements, and sample selec- 
tion. In Section 3, we will investigate the relationship 
between nuclear LINERs and extended LINERs and de- 
rive an average emission line surface brightness profile. In 
Section 4, we will show the line ratio gradient. In Section 

5, we will present the line width differences. In Section 

6, we discuss the viability of each ionization mechanism 
in explaining the line ratio gradient and line width dif- 
ferences. We conclude in Section 7. 

Throughout this paper, we use a flat ACDM cosmol- 
ogy with Q m — 0.3, and a Hubble constant of Hq = 
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Fig. 1. — Line ratio diagnostic diagram for SDSS galaxies (gray 
scale) at 0.09 < z < 0.1 with all four emission lines detected 
at more than 3cr significance, and the passive red galaxies (black 
points) among them. Th e lat ter is selected according to the cri- 
teria desc ribed in Sect ion 12.41 The curve s represent demarcations 
defined bv lKewlev eTaTI 1 1200 ll) (solid) and lKauffmann et alj 120031) 
(dashed) . This illustrates that red galaxies have fairly uniform line 
ratios which puts most of them in the LINER-like galaxy category. 



75ft.75kms 1 Mpc 1 with /175 = 1. This Hubble constant 
is chosen to make the compar ison easier with the data 
presented bv lHo et all (|1997aD . All the magnitude used 
are in the AB system. 

2. DATA AND MEASUREMENTS 

To investigate the luminosity and line ratio gradient, 
we compare line luminosity and line ratio measurements 
between different physical apertures for the same popu- 
lation of galaxies. The Palomar survey (IHo et al.l 119951) 
provides us the best sample for nuclear aperture mea- 
surements. The SDSS survey can provide consecutively 
larger aperture measurements if we select identical sam- 
ples at consecutively higher redshifts. 

2.1. Data 

In the Palomar survey, nuclear spectrum were taken 
for a sample of ~ 500 local galaxies, selected from 
the Revise d Shapley-Ames Cata l og of Bright Galax- 
ies (RSA; iSandage fc Tammannl Il981f ) and the Sec- 
ond Reference Catalogue of Bright Galaxies (RC2; 
Ide Vaucouleurs et al.1 119761 ) with the criteria of Bt < 
12.5 (Vega magnitudes) and S > 0. The nuclear regions 
(~ 200 pc) of these galaxies are isolated using a 2" x 4" 
aperture. The details of data reduction, stellar contin- 
uum subtractio n, and line measurements can be found in 
IHo et al.l dHH). 

We also employ data from the Sloan Digita l Sky Sur- 
vey (lYork et all MM, iStoughton et al.l l200l Data Re- 
lease Seven ( Abazaiian et al.ll2009D . Using two fiber-fed 
spectrographs on a dedicated 2.5-m telescope, SDSS has 
obtained high quality spectra for roughly half a million 
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galaxies with r < 17.77 in the wavelength range of 3800- 
9200A with a resolution of R ~ 2000. The SDSS fibers 
have a fixed aperture of 3" diameter, which corresponds 
to different physical scales at different distances. 

The SDSS spectroscopic data used here have been 
reduced through the Princeton spectroscopic reduction 
pipeline, which produces the flux- and wavelength- 
calibrated spectraQ The redshift catalog of galaxies used 
is from the NYU Value Added Galaxy Catalog (DR7) 
(jBlanton et all 120051). K -corrections for SDSS were de- 
rived using iBlanton fc Roweil (|2007l) 's kcorrect package 
(v4_2). 

2.2. Emission line measurements 

For emission line measurements in the Palo mar sam- 
ple, we adopt the tabulated values provided bv lHo et all 
(|1997afl . 

For the SDSS sample, we measured the emission lines 
in the spectra after a careful subtraction of the stellar 
continua. T he code us e d is a n updated version of the 
code used bv lYan et all ([2006D . The major updates are: 

1. We apply an additional flux calibration to all of 
SPSS spect ra to fix small scale calibration residuals 
(|Yanl 120111 ) . This is critically important for our 
results. We describe this correction in more detail 
below. 

2. The absolute flux is calibrated for each spectrum 
by matching the synthetic r-band magnitude with 
the r-band fiber magnitude. The spectra are also 
corrected for Galactic extinction. 

3. The stellar continuum is modelled as a non- 
negative linear combination of 7 templates. The 
templates are seven simple stellar population mod- 
els with solar metallicity, with ages of 0.125, 
0.25, 0.5, 1, 2, 7, and 1 3 Gyrs, made using the 
IBruzual fc Chariot! ()2003[ ) stellar population mod- 
els. 

As in lYan et all (|2006l ) . the emisson line flux is 
measured from direct flux-summing in the continuum- 
subtracted spectra. The line windows and sidebands are 
unchanged. 

We discovered that the flux calibration produced by 
the standard SDSS pipeline has percent-level small-scale 
residuals that can significantly impact the measurement 
of emission line flux when the equivalent width of the line 
is low (a few Angstroms; see lYanll2011l Fig. 1 for an ex- 
ample of the impact). For example, for an emission line 
EW of lA measured in a 20Awindow, if the throughput 
differs by 1% between the central window and the side- 
bands where the continuum level is measured, the line 
flux measurement will be off by 20%. This will intro- 
duce systematic offsets in line flux and line ratios as a 
function of redshift, significantly hampering our investi- 
gation. Therefore, we need a much more accurate small 
scale flux calibration than what the standard pipeline 
p roduces. 

lYanl (|2011f ) solved this problem by comparing stacked 
red-sequence spectra between small redshift intervals to 

1 http://spectro.princeton.edu/ 

2 http:/ /sdss. physics. nyu.edu/vagc/ 



statistically determine the relative throughput as a func- 
tion of wavelength, and achieved an flux calibration ac- 
curracy of 0.1% on scales of a few hundred Angstrom. 
We applied this small scale flux calibration to the spec- 
tra. This calibration is ess ential for t he result presented 
in this paper (see Fig. 6 in IYanll201lD . 

In this paper, we also make use of the line width mea- 
surements for the SDSS sample. The line widths are 
measured by fitting Gaussians to each emission line. Dif- 
ferent emission lines are allowed to have different widths, 
except that the two [S II] lines (6716A and 673lA) and 
the two [N II] lines (6548A and 6584A) are both forced 
to have the same width. The instrumental resolution of 
SDSS varies with wavelength, the position of the fiber 
on the focal plane, and the spectrograph. This varying 
resolution as a function of wavelength is given for each 
individual spectrum by the Princeton pipeline. We sub- 
tracted quadratically the instrumental broadening from 
the measured line width to derive the intrinsic width of 
each emission line for each galaxy. Our quoted uncer- 
tainty of the line width measurement is the formal un- 
certainty of the Gaussian fit. 

2.3. Photometry 

To identify the same population of galaxies at different 
rcdshifts, we use a photometric selection. We intention- 
ally avoid the use of line ratios in sample selection to 
avoid bias on the line ratio gradient estimates. 

For th e Paloma r surve y, we took the catalog pro- 
vided by IHo et al.l (|1997af ) . Photometric information is 
available from the Third Reference Catal ogue of Bright 
Galaxies fRC3 lde Vaucouleurs et al.lll991l ) but is incom- 
plete. To increase the sample size with available photom- 
etry and to put them on the same system as the SDSS 
galaxies, we re-measured photometry for those galaxies 
inside the SDSS footprint using the SDSS images. We 
employed an improve d background subtraction technique 
( Blan ton et all 120111 ) to treat these nearby large galax- 
ies properly. After proper background subtraction, mo- 
saicking, and deblcnding, we measure the photometry by 
fitting a two-dimensional Sersic profile to the deblended 
galaxy image. In the end, we derive the Galactic extinc- 
tion corrected restframe B and V magnitudes for these 
galaxies from the measured g and r magnitudes using 
the k correct software package (v4_2, IBlanton fc Roweisl 
12001 . 

For those galaxies outside the SDSS footprint and for 
certain Messier objects for which the new method does 
not yield satisfactory results, we take the photometric 
information from the RC3 catalog and convert them to 
the AB system, and then correct for Galactic extinction. 
We do not attempt to correct for internal extinction for 
these galaxies as such measurements are not available for 
the SDSS sample. 

For the higher redshift SDSS spectroscopic sample we 
derived the B and V magnitudes from the SDSS magni- 
tudes using the kcorrect package mentioned above. 

2.4. Sample definition 

Figure[5]shows the color-magnitude distribution for the 
Palomar sample overlaid on a sample of SDSS galaxies 
between 0.09 < z < 0.1. The two samples have consistent 
color-magnitude distributions. We select only the red- 
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Fig. 2. — Color-magnitude distribution of the Palomar sample 
(dark crosses) and galaxies in SDSS with 0.09 < z < 0.1 (contour 
and gray points). The lines indicate our color and magnitude cuts. 



sequence galaxies in both samples using two stringent 
color cuts denned by 

(B — V) > -0.016(M V -51og/i 75 ) +0.415 (1) 
(B-V)< -0.016(M V - 5 log h 75 ) + 0.475 (2) 

These cuts are chosen to reduce contamination from 
dusty star-forming galaxies. We limit to galaxies brighter 
than —20.4 in My — 5 log /175 to match the magnitude 
limit of the SDSS survey at z ~ 0.10. 

There are 86 red galaxies in the Palomar survey satis- 
fying these two cuts. Based on the morphological type 
given by the RC3 catalog, there are 30 ellipticals, 30 
lenticulars, 15 early spirals(S0/a, Sa, Sab), 9 late spirals 
(Sb,Sbc,Sc) and 2 irreg ular galax i es. Wi th the classifi- 
cation scheme given in IHo et al.l (j!997aft . there are 29 
LINER nuclei, 19 transition objects, 11 Scyferts and 3 
HII regions. The remaining 24 objects have no line emis- 
sion detectable in their nuclei. 

In SDSS, we select a comparison sample with 0.01 < 
z < 0.1 using the same color and absolute magnitude 
cuts. In most of our analysis, we bin the SDSS sample 
into 9 redshift bins with Az = 0.01. 

Despite our selection on color, red galaxies can have 
sufficient star formation to contribute to the line emission 
in our apertures, especially for the more distant galax- 
ies. The morphological distribution of the Palomar red 
galaxy sample suggests that this is occuring, given the 
presence of late-type spirals. We would like to exclude 
star-forming spectra in our analysis, since we want to un- 
derstand the origin of line emission not associated with 
star formation. Therefore, in the Palomar sample, we ex- 
clude galaxies with Hubble types later than SO and those 
spectroscopically classified as HII nuclei. This removes 
31% of the Palomar red sample. The remaining sample 
includes 19 LINERs, 13 transition objects, 6 Seyferts, 



Fig. 3. — D„(4000) vs. log [N II] /Ha for red sequence galaxies 
in SDSS with 0.09 < z < 0.1 and M v < -20 .4. We show only the 
brightest 50% of the sample in Ha luminosity. We show our chosen 
threshold as the solid horizontal line; it is set at the 30-th percentile 
in the D n (4000) distribution of the whole sample. Those galaxies 
with low D n (4000) and low [N II] /Ha probably have significant 
contribution by young massive stars in the production of their line 
emission. 

and 21 quiescent galaxies. 

To achieve a similar selection in the SDSS sample, we 
exclude galaxies with any star formation using a strin- 
gent cut on D„(4000) (|Balogh et al.lll999D . This quantity 
is a proxy for the light weighted mean stellar age, and 
is thus sensitive to small levels of star formation. Mea- 
sured over two 100A windows separated by 50A, it is less 
affected by dust reddening than rest frame colors, and 
more robustly measured than the H5a equivalent width. 
Figure [3] shows the D„(4000) vs. log [N II] /Ha for those 
galaxies in the SDSS sample with 0.09 < z < 0.1. Galax- 
ies that have small Z?„(4000) also tend to have lower 
[N II] /Ha, suggesting that star formation could be con- 
tributing significantly among these. In the SDSS sample, 
we remove the 30% galaxies with the lowest -D n (4000). 
This selection by D„(4000) rank is done separately in 
each redshift bin to take into account any potential aper- 
ture effects and redshift evolution. We choose a cut on 
-D„(4000) rather than one based on line ratios to avoid 
biasing the comparison of emission line properties. In 
the rest of this paper, we refer to samples with possible 
star-forming galaxies removed as the Palomar red sample 
and the SDSS red sample. 

To summarize, from the Palomar survey and SDSS, we 
identified a volume-limited sample of passive red galaxies 
without any star formation at < z < 0.1. 

3. SPATIAL DISTRIBUTION OF LINE EMISSION 

In this section, we investigate the spatial distribution of 
line emission. The spatial distribution alone does not dis- 
tinguish between different ionization mechanisms, but it 
is essential for the interpretation of other measurements, 
such as line ratio gradients. 
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We do this in a two-step process. First, we compare the 
nuclear aperture measurements from the Palomar survey 
with the large aperture measurements from SDSS at z ~ 
0.1 to establish the relation between nuclear LINERs in 
Palomar and the LINER-like galaxies in SDSS. Then, 
we utilize all apertures available to us from < z < 0.1 
to measure the average emission line surface brightness 
profile among passive red galaxies. 

3.1. Nuclear Emission vs. Extended Emission 

In this section, we will compare the emission line lu- 
minosity distributions between two identically-selected, 
volume-limited samples, but for which the line luminosi- 
ties are measured from different physical apertures. The 
difference in their luminosity distributions demonstrates 
that the emission measured in the larger aperture has to 
be spatially extended. 

In the left panel of Figure 0] we compare L(Ha) and 
[N II] /Ha between the Palomar red sample and the SDSS 
red sample at z ~ 0.1. Not all galaxies in either sample 
have Ha detected (64.4% of the Palomar red sample and 
52.0% in the SDSS red sample at z ~ 0.1 have Ha de- 
tection). Therefore, we only compare the brightest half 
of each volume-limited sample in Ha luminosity. 

The SDSS red sample shows much brighter Ha lumi- 
nosities and slightly lower [N II] /Ha ratios than the Palo- 
mar sample. Since both samples are volume-limited to 
the same magnitude cut, we adopted the same sample se- 
lection, and any evolution effect over a redshift difference 
of 0.1 should be tiny, this difference in Ha luminosities 
must be caused mostly by the difference in the physi- 
cal aperture size between the Palomar survey and SDSS. 
The Palomar sample reflects the nuclear properties of red 
galaxies while the SDSS sample reflects their integrated 
properties on much larger scales. 

The brightest 25th percentile in Ha luminosity for 
the z ~ 0.1 SDSS sample (including non-detections) is 
5.82 x 10 39 erg s _1 hy 5 2 , nearly 7 times larger than that 
in the Palomar sample (0.85 x 10 39 erg s" 1 hy 5 2 ). In fact, 
even the median Ha emitter in the SDSS red sample is 
brighter than the majority of nuclear Ha emitters in the 
Palomar red sample. Nuclear emission in red galaxies is 
therefore only rarely as luminous as found in the SDSS 
galaxies. Furthermore, we expect no strong evolution 
in AGN activity between z ~ 0.1 and z ~ 0. Thus, in 
the SDSS galaxies, a substantial contribution to emission 
must come from outside the nucleus, and therefore the 
Ha emission observed by SDSS in these z ~ 0.1 passive 
red galaxies has to be spatially extended. 

One might wonder if the spatially extended emission 
found in large aperture measurements is due to low-level 
star formation in these galaxies. We can simulate the 
expected [N II] /Ha ratio and L(Ha) by adding a typical 
star-forming emission-line spectrum to a typical nuclear 
spectrum in the Palomar sample. We use the median 
[N II] /Ha and median Ha luminosity in the brighter half 
(in Ha) of the Palomar red sample to represent a typical 
nucleus. For star formation, we adopt an [N II] /Ha ratio 
of 0.45, typical of a high metallicity star-forming galaxy, 
which yields a conservatively high [N II] /Ha ratio. The 
result is shown by the curve in the left panel of Fig. fJJ 
From the bottom end of the curve to the top, the Ha 
luminosity contributed by the star formation goes from 



to 100 times that of the typical nucleus. The curve misses 
the majority of the red galaxies at z ~ 0.1. Clearly, 
the spatially extended line emission in red galaxies we 
selected at z ~ 0.1 is not powered by star formation. 

As described above, in constructing this passively- 
evolving red galaxy sample, we have removed potential 
star-forming contaminants by removing the 30% galaxies 
with the lowest L>„(4000). The right panel of Figure[4] 
shows the effect of this procedure, where we plot the low 
D„(4000) galaxies plus those Palomar red galaxies with 
Hubble types later than SO for comparison. In the SDSS, 
those red galaxies we have removed generally have higher 
Ha luminosities and lower [N II] /Ha than those we have 
kept do. For the right panel, the curve indicates the 
track traced by adding star formation to a typical pas- 
sive red galaxy in the SDSS sample from the left panel. 
From the bottom end of the curve to the top, the Ha 
luminosity contributed by star formation goes from to 
10 times that of the median Ha luminosity in passive 
red galaxies at this redshift. The curve traces the dis- 
tribution fairly well, suggesting that the line emission in 
red galaxies with low D„(4000) have more ongoing, low- 
level star formation than those red galaxies with high 
£>„(4000). 

The right panel of Fig. 0] also demonstrates that the 
Palomar red galaxies we have removed from the sam- 
ple have fairly similar line luminosities and line ratios to 
those Palomar red galaxies we have kept (except for the 
three HII nuclei, with the lowest [N II] /Ha ratios). This 
result reflects the fact that the Palomar spectra have 
smaller physical apertures. Therefore, although the cri- 
teria we used to remove star-forming contaminants dif- 
fers slightly for the Palomar sample than for the SDSS 
sample, the Palomar sample is insensitive to these differ- 
ences. 

Now, we have shown that the line emission regions in 
passive red galaxies in SDSS have to be spatially ex- 
tended, simply because few nuclear regions in red galax- 
ies are luminous enough to explain the SDSS results. The 
next question is which galaxies host these extended Ha 
emission regions at z ~ 0. Are they the same galax- 
ies that host those nuclear emission regions? The answer 
must be "yes." Suppose they were not the same galaxies: 
then galaxies hosting these extended Ha emission regions 
would need to have undetectable line emission in their 
centers. No such galaxi es are found b y the SAURON sur- 
vey. As shown bv lSarzi et alJ (|2006l ). in a representative 
sample of 48 early-type galaxies in the nearby universe, 
all galaxies with emission lines detectable have their line 
flux peaking at the center, and the distribution is nearly 
always extended. Therefore, we conclude that most, if 
not all, red galaxies that have nuclear emission line re- 
gions also have extended line-emitting regions, and vice 
versa. The host galaxies of nuclear line emitting regions 
and those of extended line emitting regions are largely 
the same population. 

3.2. Emission line surface brightness profile 

In this section, we compare the emission line luminosity 
distributions among passive red galaxies at a series of 
rcdshifts, which translates to a series of apertures, to 
investigate the average surface brightness profile of their 
line emission. We bin the SDSS sample with 0.01 < z < 
0.1 into 9 redshift bins with a binsize of 0.01. We limit the 
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Fig. 4. — Left: logL(Ha) vs. log[Nll]/Ha for passive red galaxies in the Palomar sample (dark crosses) and an SDSS comparison 
sample at z ~ 0.1 (gray points). Only the brighter half (in Ha luminosity) of each sample is plotted, since the fainter half mostly have 
Ha undetected. The Palomar data is measured from nuclear aperture spectra while the SDSS data is measured from spectra integrated on 
much larger scales (5 kpc diameter). The curve represents the track followed by adding star formation to a typical Palomar nucleus. 
Clearly, star formation is not the cause of the luminosity difference. The much larger Ha luminosities shown by the SDSS sample suggest 
that the line emission in them is spatially extended. The right panel shows the red galaxies that we have excluded from the left panel 
because of possible star formation contamination. For the SDSS, this means those red galaxies with the lowest D n (4000) ; for Palomar, this 
means those red galaxies with Hubble types later than SO and those classified as HII nuclei by Hp et al. (1997aj). The curve represents the 
track followed by adding star formation to a typical passive red galaxy in SDSS from the left panel. For the SDSS samples in both panels, 
we also excluded a small fraction (11% of the bright half in the left panel, and 1.5% in the right panel) of galaxies with poorly measured 
[N Il]/Ha ratios (with uncertainties larger than 0.25 dex), which has no practical impact on our conclusions. 



Palomar sample to only those galaxies at D < 40Mpc, 
which corresponds to z = 0.01, for the lowest redshift 
bin. 

Figure [5] shows the Ha luminosity distribution as a 
function of redshift. In each bin, we only plot the brighter 
half of the sample in Ha. With increasing redshift, i.e., 
increasing aperture, the Ha luminosities increase. Thus, 
the Ha luminosities observed with larger apertures have 
to have significant contributions from spatially extended 
emission line regions. In each redshift bin, we sort all 
passive red galaxies by their Ha luminosity. Figure [6] 
plots the brightest 25th percentile Ha luminosity as a 
function of physical scale covered by the SDSS fiber. 
The 25th percentile is safely above the detection thresh- 
old at all redshifts. The luminosity increases with scale 
roughly as a power law with an index of 0.72, as shown 
by the power-law fit in Fig. [6] As this is the integrated 
luminosity within radius r, it indicates the average sur- 
face brightness profil e follows r~ 1 28 . T his slope is fairly 
consistent with what ISarzi et al.l (|2010f ) found in nearby 
early-type galaxies targeted by the SAURON survey (see 
their Figure 4). We also show the same measurement of 
the Palomar sample (with late-type galaxies removed) at 
the median effective radius probed by the 2" x 4" aper- 
ture in the Palomar survey, which we treat as equivalent 
to a circular aperture with 3"diameter. The 25th per- 
centile in the Palomar sample is fairly consistent with the 
power-law fit to the SDSS data points. This evidence fur- 
ther strengthens the conclusion that the nuclear and the 
extended line emitting regions exist in the same galaxy 



population. 

However, the extended line emission alone does not 
constrain the source of the ionizing radiation. Contrary 
to intuition, the extended emission could also be pro- 
duced by a central ionizing source, such as an AGN. The 
emission line brightness profile depends on many factors: 
the ionizing flux profile, the density profile, the gas filling 
factor, and how the typical size of the gas clouds change 
with radius. We leave the detailed calculations to £16.11 

4. LINE RATIO GRADIENT 

With the above technique, we can also check if 
the line ratio distribution in this population changes 
between redshifts/apertures. This check can only be 
done on those galaxies with detectable line emission. To 
ensure low uncertainty on the line ratio measurement, 
we choose only the brightest 25% in total emission 
line luminosity at each redshift and compare their 
various line ratios. To avoid bias on the line ratios, 
instead of selecting the brightest 25% in Ha luminosity, 
we base the selection on the total luminosity of the 
several strongest emission lines available in the spectra 
(Ha+[NII] A6583+[0 III] A5007+[S II] AA6716,6731). 
This combination is a better proxy of the total line 
emission output than L(Ha). 

Figure [7] shows the full distribution of several line ra- 
tios as a function of aperture radius. We also plot those 
galaxies in the Palomar red sample with D < 40 Mpc to 
probe the smallest scales. As for the SDSS sample, we 
only select the brightest 25% galaxies in total observed 
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Fig. 5. — Ha luminosity distributions of passive red galaxies as a 
function of redshifts. The first bin at z < 0.01 is from the Palomar 
sample and the rest are from the SDSS sample. Only the brighter 
half of the sample in each redshift bin is plotted. The brightest 10 
percent of galaxies in Ha in each bin are plotted as points. The gray 
scales indicate the ranges populated by different percentiles in each 
bin: brightest 10-20th, 20-30th, 30-40th, and 40-50th, from top to 
bottom, respectively. The dashed line at the bottom indicates 
the 3<x detection limit in SDSS, which is 3 times the luminosity 
corresponding to the median Ha flux error in each bin. Clearly, 
the Ha luminosity increases with increasing redshift or aperture 



emission line luminosity. 

Fig. [7] shows the distribution of line ratios in the red 
galaxies. At small aperture radii, the scatter is domi- 
nated by intrinsic variations in line ratios among galax- 
ies. At the large aperture end, the uncertainty in line 
measurements starts to dominate the scatter, leading to 
a slight increase in the width of the distribution. In most 
line ratios, the intrinsic variation has a standard devia- 
tion of approximately 0.1 dex. 

Interestingly, in some line ratios, the median of the 
distribution changes systematically with aperture radius, 
most notably in [N II] /Ha, [S II] /Ha, and [O III]/[S II]. 
The trends also extend to the Palomar sample at the 
smallest scales. The changes are so large in [S II]/Ha 
and [O III]/ [S II]ratios that most of the Palomar sample 
populates only one side of the median of SDSS sample 
even in its first bin. 

Figure [8] shows how the median line ratios change with 
aperture size, with the error bars giving the uncertain- 
ties of the median estimates. The systematic changes 
in [N II] /Ha, [S II]/Ha, and [O III]/[S II] with radius are 
significant and the Palomar sample confirms the trend 
on small scales. 

One might worry that the change in [O III] / [S II] ra- 
tio could be due to dust extinction, because these two 
lines are separated in wavelength. However, the Ha/H/? 
ratio is nearly constant with radius, indicating that the 
average dust extinction does not vary much. The values 
of the median Ha/H/3 ratio are also close to the Case B 
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Fig. 6. — The 25th percentile (counting from the brightest) Ha 
luminosity (asterisks with error bars) measured with SDSS fibers in 
a sample of non-star-forming red sequence galaxies as a function of 
the physical scale probed by the fiber. The triangle point represents 
the measurement in the Palomar sample, which is fairly consistent 
with the extrapolation of the power-law fit on small scales. The 
uncertainties of these measurements are measured with bootstrap 
resampling. The '+' signs at the bottom indicate the 3<j detection 
limits. 



prediction of 2.85 (or 3.1 if collisional excitation is in- 
cluded), indicating that the level of dust extinction is 
low. Applying extinction corrections to each galaxy in 
the sample only shifts the median [O IH]/[S II] ratio up 
by a nearly constant ~ 0.1 dex in each bin. 

The line ratios presented in Fig. [8] are cumulative mea- 
surements: they reflect the luminosity-weighted average 
line ratios within aperture radius r, rather than that in 
an annulus at radius r. We need to combine these median 
line ratios in integrated apertures with the median line 
luminosity of this sample in corresponding apertures, to 
derive the average line ratios in circular annuli at radius 
r. 

In Fig. [§] we show the median line luminosities of this 
sample in corresponding apertures for Ha, [S II], and 
[O III]. This plot is similar to Fig. [6l but uses only the 
25% brightest galaxies in total emission line luminosity. 
From this figure, it is evident that the [S II] line has a very 
different surface brightness profile from Ha and [O III] , 
consistent with the trends in integrated line ratios. 

With the integrated line ratios and integrated lumi- 
nosities, we can compute the line ratio in annuli. For 
example, the average [N II]/Ha ratio between radius 
and Tj can be computed by the following equations, 



[N H] 
Ha 



[NH] 
Ha 



%(Ha) ~ 



r<rj 



[N II] 
Ha 



Li (Ha) 



Lj(Ha) - Li (Ha) 



.( 3 ) 

where L^(Ha) and Lj(Ha) are the median luminosities 
in apertures with radius r, and rj, respectively. The 





Fig. 7. — Integrated line ratio distribution as a function of the aperture radius for passive red galaxies in SDSS (gray points) and 
the Palomar sample (red crosses). The line ratios are "integrated" in the sense that they are derived using the fluxes interior to the 
indicated aperture radii. Only the 25% brightest (in Ho) galaxies in each redshift bin are plotted to ensure small uncertainty on the line 
ratio measurement. The lines indicate the 5-, 25-, 50-, 75-, 95-percentile points in each redshift bin. The aperture of Palomar galaxies are 
computed using an effective radius of 1.5". 



[S II] /Ha and Ha/H/3 ratios in annuli are also computed 
by combining integrated line ratios with the median Ha 
luminosities. The [N II] /[S II] and [O III]/[S II] ratios 
are computed by combining integrated line ratios with 
the median [S II] luminosities. The [O III]/H/3 ratios are 
computed by combining with the median [O III] lumi- 
nosities. Calculating this between every consecutive bin 
results in large uncertainties, due to the large fractional 
error in the luminosity differences. We therefore calcu- 
late the annulus line ratios using aperture pairs: = 
[1, 3], [1,4], [2, 5], [3, 6], [4, 7], [5, 8], [6, 9], [7, 10], and[8, 10], 
where Aperture 1 is the Palomar aperture and Aperture 
10 is the aperture for SDSS galaxies at 0.09 < z < 0.1. 
The results are shown in Fig. [TU] 

Although the uncertainties become much larger, the 
differences between the innermost bin and the outer bins 
remain robust in [N II] /Ha, [S II] /Ha, and [O III]/[S II]. 
The second bin in each panel always shows a very large 
uncertainty. This is due to the dramatic change in line 
ratios between the first few bins and their larger uncer- 
tainties. In Table [1] we list the line ratios of the inner 
most bin and the median line ratios of the outer bins 
(700 < r < 2500 pc). 

The line ratio gradients we observe have been in prin- 
ciple detectable in past l ong-slit spectroscopic surveys of 
early-type galaxies (e.g.. iPhillips et al.|[T9 86; Ki mlll989l : 
IZeilinger et al]|1996t iHo et al.lll997at iCaon et al.ll2000D . 
However, these authors either did not have data with suf- 
ficient quality or did not lo ok at the line ratio gr adient 
at all. The only exception is IZeilinger et a l. (1996), who 



TABLE 1 

Median line ratios within 300 pc radius and outside 700 pc 

Location log[N Il]/Ha log[S Il]/H« log [O HI]/ [S II] 

Inner 300 pc 0.25 ±0.03 0.13 ±0.04 -0.40 ± 0.04 

Outside 700 pc 0.01 ± 0.01 -0.07 ±0.01 -0.07 ± 0.01 



showed that the [N II] /Ha ratio decreases outwards in 
four galaxies. However, like many of these past surveys, 
their data did not have wide enough wavelength cover- 
age to cover multiple line ratios, which would be critical 
for ide ntifying the cau s e of t he line ratio gradients. Re- 
cently, Annib ali et al.l (]2010l ) investigated the line ratio 
gradients using long-slit spectra with wide wavelength 
coverage for a sample of relatively gas-rich early-type 
galaxies. They found that the [N II] /Ha ratios in most 
of them decrease with radius, consistent with our result. 
However, they did not look for gradients in other line ra- 
tios, except for [O III]/H/3, which is very often too noisy 
to support firm conclusions. 

We will discuss in $6] what changes in physical condi- 
tions are required to produce such changes in line ratios. 

5. CLUES FROM LINE WIDTHS 

We initally thought that the line ratio gradient would 
produce different line widths in different emission lines, 
due to varying kinematics of the line emitting clouds 
across each galaxy. If this were so, it could provide a 
complementary constraint on the radial distribution of 
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Fig. 8. — Median integrated line ratios as a function of aperture radius for the 25% brightest line emitter s among the red g alaxies in 
SDSS (> 300 pc) and the Palomar sample (< 300 pc). For the median, we use the biweight center estimator (Beers et al. 1990), which is 
more robust than the median for small samples. The uncertainty is given by the jackknife error of the biweight estimator. All panels cover 
0.5 dex in y-axis for the ease of comparison. The horizontal line in Hct/H/3 panel marks the Case B ratio of 2.85. 
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Fig. 9. — The median integrated Ha (asterisk), [S II] (triangle), 
and [O III] (square) luminosities as a function of aperture among 
the 25% passive red galaxies (in each bin) with the brightest total 
emission line luminosity. The smallest scale is probed by the Palo- 
mar sample and the larger scales are probed by the SDSS sample. 
The [S II] and [O III] points are slightly shifted in the horizontal 
direction for clarity. The uncertainties of these measurements are 
measured with bootstrap resampling. 



the emission. The data do indeed show different line 
widths for different lines. However, we have concluded 
that this variation is probably not due to line ratio and 
kinematic gradients, primarily because the line width dif- 
ferences are not a function of aperture size. We describe 
our investigation of line width differences in this section. 

Because the width measurement is noisy on SDSS spec- 
tra, we need a control sample to demonstrate that our 
line width measurement is robust. The star-forming 
galaxies provides such a comparison. In a pure star- 
forming galaxies, the line emission originates mostly from 
HII regions photoionized by O and B stars. The line 
width in integrated spectra reflects the rotation velocity 
of the galaxy. In the absence of strong metallicity gradi- 
ent, we would observe approximately the same line ratios 
in all HII regions. In this case, all lines should display 
the same line width. 

We choose a star forming g alaxy sample using the line 
ratio criteria described by IKewlev et al.l (|2006f) . Basi- 
cally, the galaxies in this sample are selected to fall in 
the star-forming branch on all three diagnostic diagrams 
([O III]/H/3 vs. [N II]/Ha, [S II]/Ha, and [O I]/Ha). We 
compare this star-forming galaxy sample with the sam- 
ple we used for deriving the line ratio profile, namely the 
top 25% passive red galaxies that have the highest total 
emission line luminosities in each redshift bin. We com- 
bine all the redshift bins together. In addition, to ensure 
good measurements on the line width ratio, we require 
the uncertainty of the line width ratio to be smaller than 
0.1 dex. 

In Figure 111! we plot the distributions of width ratio 
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Fig. 10. — Average line ratios in annuli as a function of radius for passive red galaxies. This is derived by combining the integrated 
line ratios in Fig. [8] with the integrated luminosities in Fig. [9] through equations similar to Eqn.[3] The horizontal bar of each data point 
indicates the range of radius covered by the annulus. Note the y-axis in each panel covers a larger range (0.7 dex) than Fig. [8] 



in four line pairs for galaxies in the star-forming sam- 
ple (thin-lined histograms) and in the passive red galaxy 
sample (thick-lined histograms). In all line pairs, the 
distribution of star-forming galaxies always show a fairly 
symmetric and narrow distribution peaking around zero 
in logarithmic space, meaning all the lines have roughly 
the same width. This proves that our line width mea- 
surement is robust. It also indicates that line ratios in 
star-forming HII regions do not correlate strongly with 
the velocity of the HII regions. The [O III]-to-[S II] pair 
may be an exception; for star-forming galaxies this width 
ratio has a wider distribution than the other line pairs. 
This broad distribution probably reflects intrinsic varia- 
tion in the population. It might be caused by variation in 
the disturbed component of the diffuse ionized medium 
in star- formi ng galaxies, which produces strong and wide 
[O III] lines (jWang et all 119971) . We leave this question 
for future investigations. 

However, for line-emitting passive red galaxies, the dis- 
tributions in [N II]-to-[S II], [O III]-to-[S II], [S II]-to-Ha 
width ratios do not peak around zero in log space. Their 
offsets from zero are highly significant. On average, the 
[N II] and Ha lines in them are wider than [S II] lines 
by 8% and the [O III] lines are wider than [S II] by 16% 
(Table [3]). This means the [O III]/[S II] ratios in the line 
wings are higher than that in the line center, and the 
higher velocity clouds which contribute to the wings must 
have higher [O III] / [S II] ratios. The same must be true 
for [N II] /[S II] and Ha/[S II] ratios. This indicates that 
in these galaxies, line emitting clouds do not have uni- 
form line ratios and the line ratio must correlate with the 



velocity of the clouds. 

In the previous section, we found the line ratios have a 
systematic variation with radius. Could these line ratio 
variations produce the line width differences? 

Given the fact that [O III] / [S II] ratio increases with 
radius, to obatin a wider [O III] line than [S II], the line- 
of-sight velocity broadening must also increase with ra- 
dius. The broadening could either come from random 
motions or ordered rotation. As sh own by previous long 
slit (IPhillips et alJfl^lK^lrMllZeihnger et aLlUflflfift 
and IFU (|Sarzi et ah 2006) observations, in most early- 
type galaxies, the kinematics of the gas is largely con- 
sistent with disk-like rotation with rotation velocity in- 
creasing outwards (though see iHeckman et alJ 119891 for 
some exceptions). This is consistent with the require- 
ment here. The [N II]-to-[S II] and [S Il]-to-Ho width 
ratios could also be consistent with their respective line 
ratio gradients (if the decrease in [N II] / [S II] towards the 
center is real, as indicated by the Palomar data point). 
However, the [N II] show the same width as Ha, incon- 
sistent with the expectation. 

On the other hand, if the width difference is indeed 
caused by the line ratio gradient and the rotation of the 
gas disk, the line width difference should gets smaller in 
larger apertures, since we expect the flux to be increas- 
ingly dominated by the outskirts where the line ratio 
profile gets flat. Figure [T2l shows the line width ratio dis- 
tributions as a function of aperture radius. Apparently, 
the average line width ratios are roughly constant, inde- 
pendent of the aperture size. This is inconsistent with 
the expectation of the model. 



12 



Yan et al. 



0.20 



0.15 



c 

o 



a 0.10 
u 

tu 



0.05 



0.00 



-0.4 



Passive Red Galaxies 
Star-forming Galaxies 

J 




-0.2 0.0 0.2 

L °g ( CT [N II]/ CT [S II]) 



0.10 



0.08 



0.06 



0.04 



0.02 



0.00 



Passive Red Galaxies 
Star-forming Galaxies 




-0.2 0.0 0.2 
L °g ( CT [o ni]A[s ii]) 



0.4 



0.25 



0.20 



o 0.15 

o 
co 

0.10 



0.05 



0.00 



-0.4 



Passive Red Galaxies 
Star-forming Galaxies 




-0.2 0.0 0.2 

Log (ff [NII] /CT Ha ) 



Passive Red Galaxies 
Star-forming Galaxies 




-0.2 0.0 0.2 
Log (ff [s n] /a Htt ) 



0.4 



Fig. 11. — Panel (a): Distributions of line width ratio between [Nil] and [S II] for star-forming galaxies (thin-lined histogram ) and 
passive red galaxies (thick-lined histograms). [N II] in passive red galaxies are systematically wider than [S II] in both samples. Panel (b), 
(c), (d) show the same comparisons in [O IIl]-to-[S II] width ratio, [N Il]-to-Ha width ratio, and [S Il]-to-H« width ratio, respectively. 



We demonstrate this inconsistency quantitatively by 
simulating the expected line width ratios using a toy 
model of a rotating gas disk in a spherically-symmetric 
galaxy with stellar mass of 7.5 x 10 10 Mq and an effec- 
tive radius of 4.8 kpc, the median values among our 25% 
passive red galaxy sample with the brightest total line 
luminosity. The rotation curve of the gas disk is set by 
the stellar density profile, which is assumed to be a 7- 
model described by iDehnenl (| 1993T ) with 7 = 1.5. This 
stellar density profile gives a stellar surface brightness 
profile closely resembling the de Vaucouleurs' i? 1 / 4 pro- 
file. The integrated luminosity profile in [O III] is fixed 
to be a power law with index of 0.77, measured by fitting 
the data points in Fig. |H1 The inclination of the disk is 



set at 60°. We assume that the line ratios at each point 
in the disk are solely dependent on radius. We model the 
logarithm of the line ratio profile as a broken power-law 
of the form, 

[O III] = J A(r/r )^ - 0.5 , r < r , , 

B [S II] \ (A - 0.4)(r/r )^ - 0.1 ,r > r y ' 

Based on the trend seen in Fig. [101 we fix the model to 
have log [O III] / [S II] equal to -0.5 (an arbitrary choice) 
at r = and asymptote to -0.1 as r — s> 00. We fit the 
model to the integrated line ratios rather than the dif- 
ferentiated line ratios, since the former have independent 
uncertainties. 

The velocity dispersion in the disk is assumed to be a 
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Median line width ratios 
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Fig. 12. — Top panel: the [O IIl]-to-[S II] width ratio distribution 
as a function of aperture radius for the 25% passive red galaxies 
with the highest total line luminosity. The dark points with error 
bars indicate the median width ratio in each redshift bin. Bottom 
panel: same plot for the width ratio between [N II] and [S II] . 



TABLE 2 

Line-ratio Profile Model Parameters 



A 


r (kpc) 


7i 


72 


1.86 


345.9 


3.69 


-2.44 


1.56 


309.2 


3.82 


-2.04 


1.21 


402.3 


2.29 


-1.98 


0.982 


565.1 


1.30 


-2.30 


0.868 


364.7 


1.58 


-1.46 



constant everywhere and equal to 50kms _1 . Many gas 
kinematic studies have shown that the velocity disper- 
sion is likely to increase towards the center. Here we as- 
sume the extreme case of flat dispersion, since using an 
increasing velocity dispersion towards the center would 
erase the line width differences. 

We employ a Markov-Chain Monte Carlo technique to 
find a large number of models that best fit the data. 
Among the 10 data points, we ignored the second bin 
(r ~ 500 pc) in the fitting, as including it makes the 
fit difficult. We chose five typical but different models 
to illustrate the trend expected in line width differences 
as aperture radius increases. The model parameters are 
given in Table. [2] 

The left panel of Fig. [T3] shows these models give rea- 
sonably good fits to the integrated line ratio profiles. The 
right panel of Fig. [T3l shows the line width ratio as a func- 
tion of aperture radius for these models. The solid curves 
and the dashed curves show the result for two different 
methods for measuring the line width. Because circular 
rotation makes more boxy profiles than Gaussian, the 
width measured from Gaussian fitting is different from 
that measured from FWHM. No matter how the width is 
measured, the width ratio between [O III] and [S II] pro- 



Line pair 


Star-forming galaxies 


Old red galaxies 


CT [N II]/ CT [S II] 
CT [0 III]/ CT [S II] 
a [S II] Mia 
CT [N Il]/ CT Ha 


0.991 ± 0.0004 
1.006 ±0.001 
1.033 ±0.0004 
1.022 ±0.0003 


1.079 ±0.004 
1.164 ±0.006 
0.933 ± 0.004 
1.009 ±0.003 



duced by all circular rotation models decreases strongly 
towards large apertures, inconsistent with what we ob- 
serve. The decrease in width ratio towards large aperture 
is expected in the model since the line luminosity is in- 
creasingly dominated by flux from large radius and the 
line ratio profile flattens outward. Therefore, the line 
width ratio is probably not caused by the line ratio gra- 
dient and circular rotations. So, what is the real cause 
of the line width difference? 

The line width is mainly contributed by two compo- 
nents: the thermal broadening and the bulk motion of 
the clouds. For gas at T = 10 4 K, the thermal broaden- 
ing is approximately 15km/s for Ha. For the red galaxies 
in our sample, the lines are very wide, with Gaussian sig- 
mas ranging between lOOkm/s and 300km/s and a me- 
dian of 176km/s in [Nil] lines. Therefore, the thermal 
broadening is a minor contributor. The width is likely 
dominated by bulk motion broadening or turbulence in 
the clouds. 

To produce a width difference between different lines, 
we have to have clouds with different line ratios and dif- 
ferent line widths, and the line ratio has to correlate with 
the line width. For example, for [O III] and [S II] to have 
different widths, we need a population of clouds with high 
[O III] / [S II] flux ratio and a population of clouds with 
low [O III]/[S II]. To make a wider [O III] than [S II], 
those high [O HI]/[S II] clouds need to produce a wider 
width than those low [O I1I]/[S II] clouds. Basically, to 
produce a constant width ratio with radius, the width 
ratio has to be approximately the same everywhere in 
the galaxy. Thus it requires at least two components of 
the ISM throughout the galaxy that have different line 
ratios and different line broadening. Because the total 
flux ratio has a gradient with radius, the two compo- 
nents need to have approximately synchronous line ra- 
tio gradient behaviours. One component is kinemati- 
cally more disturbed than the other and thus produces a 
wider line width. Our observed width ratios indicate that 
the more disturbed component has higher [O III] / [S II] 
and [N II] /[S II] ratios, lower [S II] /Ha ratio, and similar 
[N If |/Ha ratio to the more quiescent component. With- 
out knowing the ionization mechanism for the gas, the 
origin of these multiple components and the reason for 
their line ratio difference is difficult to analyze. 

We check if the line width difference changes with the 
line strength. Figure [14] shows the width ratio distri- 
bution between emission lines as a function of Ha EW. 
We only use the brightest 25% galaxies in total emission 
line luminosity and excluded among them those with un- 
certainty on line width ratio greater than 0.2 dex. The 
median [O III]-to-[S II] and [N II]-to-[S II] width ratios al- 
ways stay above 1. The [O III]-to-[S II] width ratio seems 
to decrease slightly with increasing Ha EW. 
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Fig. 13. — Left panel: [O III]/[S II] line flux ratio as a function of aperture radius for the data (points with error bars) and a few simple 
models. Right: The [O IIl]-to-[S II] width ratio as a function of aperture radius for the data and the models. The two sets of curves indicate 
line widths measured by different methods, with solid lines indicate FWHM ratios and dashed lines indicate ratios from Gaussian fits. The 
models assume thin gas disk in circular rotation, with flat velocity dispersion. 
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Fig. 15. — Top panel: the width ratio distribution between [O III] 
and [S II] as a function of luminosity for line-emitting red galaxies. 
Bottom panel: same plot for the width ratio between [N II] and 
[S II]. 



Figure [15] show the width radio distribution for the 
line-emitting red galaxies as a function of absolute lu- 
minosity. The same sample is used as in Fig. [14] The 
[N II]-to-[S II] width ratio stays flat as a function of lumi- 
nosity but the [O III]-to-[S II] width ratio declines slowly 
towards fainter galaxies. The important point is that 
they always show significant offset from zero in log space, 
suggesting that the reason causing the width difference 



is universal in these galaxies. 

6. DISCUSSION 

In this section, we first discuss what physical fac- 
tors determine the emission line surface brightness pro- 
file, and show that the profile alone does not provide a 
discriminator between different ionization mechanisms. 
Then, we discuss the viability of the different ionization 
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mechanisms in light of the observational results we pre- 
sented above. 

6.1. Surface Brightness Profile 

Except shock heating and conductive heating by the 
hot gas, all other major ionization mechanisms proposed 
involve photoionization. In this section, we consider a 
generic photoionization model and examine which of its 
parameters determine the emission line surface bright- 
ness profile. 

We assume that the ISM is filled with hot, ionized gas. 
Embedded in it are neutral dense clouds. Each cloud is 
optically-thick to the ionizing radiation. In photoioniza- 
tion equilibrium, for each cloud the total emission line 
luminosity has to be equal to the total photoionizing lu- 
minosity it receives, which is the incoming flux times the 
projected area of the cloud. Therefore, the luminosity 
density profile depends on the photoionizing flux profile 
(as a function of radius) and the total projected cloud 
area per unit volume. For example, for a volume filling 
factor of /„, assuming the clouds have an average vol- 
ume of (V) and an average projected area of (A), the 
luminosity density of line emission would be 



{A) 



(5) 



Here, F(r) is the ionizing flux profile. The second term 
on the right hand side gives the number density of clouds. 
Multiplying it with the average projected area yields 
the total projected cloud area per unit volume. To ob- 
tain the final surface brightness profile, we also need to 
convolve the luminosity density profile with the spatial 
distribution of the clouds. If the clouds all reside in a 
disk with constant thickness, then the surface brightness 
scales with radius in the same way as the luminosity den- 
sity. However, if the thickness of the disk increases with 
radius, like the Milky Way gas disk, then the surface 
brightness profile would be much shallower. Assuming 
the scale height of the disk is H(r), the surface bright- 
ness of the line emission would be 



E(r) = F(r)f g (r)H(r) 



(V) 



(6) 



The typical cloud area and volume could also change with 
distance from the center. Because the gas density (n) 
depends on distance from the galaxy center, if the mass 
distribution of the clouds is independent of the distance, 
then the typical (A)/(V) will scale as n 1 / 3 . Therefore, 
the constraint from the surface brightness profile is 



£(r) = F(r)f g (r)H(r)n 1/3 oc r" 



1.28 



We do not have enough information about the geom- 
etry of the cloud distribution and how the typical cloud 
sizes change with radius to constrain the photoionizing 
flux profile. Therefore, the extended line emission only 
provides a partial constraint on the source of the ioniza- 
tion. We need more information from other methods. 

Many ionization mechanisms have been proposed to 
explain the observed emission line ratios in these red 
galaxies, which mostly have LINER-like line ratios. In 
the following sections we consider these mechanisms, di- 
viding them generically into three categories: a central 
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Fig. 16. — Line ratio diagnostic diagrams for the center (triangle) 
and outskirts (square) of passive red galaxies, overlaid on the MAP- 
PINGS III models of a classical dust-free AGN photoionization 
model. The models assume a metallicity of 2Zq and a hydrogen 
density of 1000 cm -3 . The solid lines indicate constant ionization 
parameters (logf/ = —4.0, —3.6, —3.3, —3.0) and the dashed lines 
indicate constant spectral indices (a = —2.0, —1.7, —1.4, —1.2). 

photoionizing source, distributed photoionizing sources, 
and shocks. 

6.2. Photoionization by an accreting SMBH 

An accreting supermassive black hole will emit X-rays 
and extreme UV radiation that photoionizes surrounding 
gas clouds and produces line emission. In this section we 
will examine the predictions of this model for the line 
ratio gradients and line width differences. 

6.2.1. Line ratio gradients for a SMBH 

First, we will demonstrate that this model cannot 
explain the observed line ratio gradients. To do so, 
we us e models calculated with t he MAPPINGS III 
codes (iDopita fc Sutherland! 119961 : iGroves et all l2004at 
I Allen et al.l I2008T ) . Figure [TBI shows two line ratio diag- 
nostic diagra ms with the g rids re presenting the models 
presented bv IGroves et al.l (|2004al ). 

This standard photoionization model (with no dust or 
radiation pressure) assumes a metallicity of 2Z and a 
hydrogen density of 1000 cm -3 . With the given range of 
parameters, this model cannot produce the [S II] /Ha and 
[N II] /Ha ratios observed in the center of these galax- 
ies. In fact, none o f the dust-free classical models in 
IGroves et al.l (|2004af ) can produce the central [S II] /Ha 
ratio: they are all too low. It may be possible to fit 
the [S II]/Ha ratios by tweaking the N/S abundance ra- 
tio, or adopting the dus ty, radiation-pressu re dominated 
photoionization model (jDopita et al.ll2002T ). However, it 
is unlikely that the latter model is applicable to the low 
intensity radiation fields in LINERs. 

Since we are not yet certain about the ionization mech- 
anism, it is premature to constrain the exact physical pa- 
rameters using these measurements. However, the mod- 
els should provide a guide as to the direction of change 
in the line ratios. Thus, we only use these models to 
investigate what the line ratio gradients tell us about 
the change in the physical parameters, not their precise 
values. 

The line ratios are primarily determined by four pa- 
rameters: gas density (n), metallicity (Z), ionization pa- 
rameter (logC/jl, and the spectral index (a). From the 

3 U is the dimensionless ratio of the ionizing photon flux density 
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Fig. 17. — [O III] / [S II] vs. [N Il]/Ha for the classical dust-free AGN photoionization models (grids) and the measurements at the center 
(triangle) and outskirts (square) of line-emitting red galaxies. We explore the dependence of line ratios on four parameters: ionization 
parameter (log U), spectral index (a), gas density (n), and metallicity (Z). In each panel, we fix two parameters and vary the other two 
to see how the line ratios depend on each parameter. 



four strong lines, [S II], [N II], Ha, [O III], we have only 
three line ratios at each position. Therefore, we cannot 
hope to determine the trend in all of these parameters 
and have to start by keeping some parameters fixed. In 
the six panels of Fig. [171 we fix two parameters at a time 
and look at the line ratio dependence on the other two 
parameters, to determine all possible scenarios for the ob- 
served line ratio variation. We only look at [O III] / [S II] 
vs. [Nil] /Ha diagram since the [S II] /Ha data points 
are not well covered by the models. 

There are three ways that the [O III] / [S II] ratio can 
increase outward: an increase in the density (panel c), a 
decrease in the metallicity (panel d) , or an increase in the 
ionization parameter (panel a). For the first option, it 
is unphysical to expect the density to increase outwards, 
which would require a higher gas pressure at larger radius 
(since the temperature in the ionized gas is likely to be 
always near 10 4 K). 

The second option, the metallicity gradient, is also not 
a likely source for the change in [O III] / [S II] ratio. At 
fixed density and ionization parameter, it requires a fac- 
tor of ~ 3 change in metallicity with radius, which is com- 

to the electron density, U = q(H°) / (cnn). 



parable though somewha t larger than the observ ed stel- 
lar metallicity gradients (|Kuntschner et ai1l2010[ ). How- 
ever, the density is likely to decrease outwards, which 
acts as a countervailing force and would require a larger 
metallicity gradient. Indeed, since at low metallicity the 
[O III] / [S II] ratio becomes insensitive to Z, this possibil- 
ity could be ruled out. 

The third option, an increasing ionization parameter, 
is much more promising. [O III] / [S II] is quite sensitive to 
log U: a variation of more than 0.1 dex would dominate 
any other possible effect. Even before considering specific 
models, we should suspect that the ionization parameter 
in these objects increases outwards. 

The ionization parameter is defined as the ratio be- 
tween ionizing flux and gas density. Therefore, we now 
examine how gas density changes with radius. X-ray ob- 
servations of giant ellipticals have shown that the hot gas 
density follows the square root of the stellar density pro- 
file, n e oc (hj" 2 (|Mathews fc Brighentil2003l and references 
there in). This means that the gas density profile falls 
with radius as r~ p with 0.5 < p < 1 at the center and an 
increasing p at large radii. This range of central density 
slope is c onsistent wi t h the more recent X-ray measure- 
ments bv lAllen et alj ([2006). The temperature profile of 
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the hot gas is much natter, varyi ng by at most 50% be- 
tween the center and the outskirts (|Mathews fc Brighentil 
2003). The gas clouds that generated the observed op- 
tical line emission always have temperature near 10 4 K. 
Therefore, under pressure equilibrium, with the nearly 
constant temperature profile in both the hot gas and the 
warm ionized gas, the density in the warm ionized gas 
clouds should fall with radius in roughly the same way 
as the hot gas density. 

It is important to note that this density scaling is only 
verified in giant ellipticals. In fainter early-type galax- 
ies it may not hold. Nonetheless, we expect the central 
density profile in faint ellipticals is also much shallower 
than 1/r 2 . Evidence for this comes from density profile 
measurements in the central regions (< lOOpc) of a few 
fainter early-type galaxies ( NGC 1052, NGC 39 98, NGC 
4579) from [S II] line ratios. I Walsh et all (|2008fl showed 
that the power-law indices of their gas density profiles 

are around —0.6. Therefore, we use the n e oc p~J 2 scal- 
ing as a working assumption here and in the next section. 
Our main conclusion remains the same if one switches to 
a power law density profile with n e oc r _1 or shallower. 

In the case of a central ionizing source, the flux de- 
creases as r~ 2 . Since the gas density is at most decreas- 
ing as r , the ionization parameter must decrease out- 
wards by a large amount, at least 0.5 dex. No change 
in metallicity or spectral index could conceivably make 
up for this decrease: this model inevitably predicts a 
strongly decreasing [OlII]/[SlI] ratio with radius, the 
opposite of what we observe. Therefore, unless the gas 
density profile actually falls faster than r~ 2 , the AGN 
photoionization model cannot explain the observed line 
ratio gradients. 

Meanwhile, there are also three ways for the [N II] /Ha 
ratio to decrease outward: a softening of the ionizing 
spectrum (panel a), a decrease in the density (panel b), a 
decrease in the metallicity (panel d) , or a combination of 
these. The outward decreasing density provides a natural 
solution. Though its decline with radius might be too 
slow to explain all the change in [N II] /Ha ratio. An 
additional contribution from metallicity gradient and a 
change in the spectral index might be needed as well. 
With the current modeling uncertainty, we cannot break 
the degeneracy among these possibilities. 



iDe Robertis fe Osterbrockl (|1986l ). the median of [N II 



to-[S II] width ratio is 1.11 ± 0.06, the median of [O III 
to-[S II width ratio is 1.18 ± 0.08. Our width ratios are 
o nly slightly smaller. 

Filippcn ko fc Halpernl (|1984| ) proposed the following 
picture to explain the line width differences. The ion- 
ization flux decreases outward according to the inverse- 
square law. If the density also falls as r -2 , then the ion- 
ization parameter seen by each cloud will be the same. 
If all clouds are optically-thick to the ionizing radiation, 
then they will all have the same ionization structure and 
produce the same set of emission lines. The relative line 
ratio will vary according to the density of each cloud. 
Lines with a high critical density will be mainly con- 
tributed by a high density clouds, which are closer to the 
nucleus and have higher velocities. 

In this scenario, the [N II] /[S II] and [O III]/[S II] flux 
ratios should increase towards the center, the opposite 
of what we observe for LINERs. In addition, for most 
red galaxies the line emission is spatially extended, with 
an average surface brightness profile falling as r^ 1 - 28 . 
Thus, the line luminosity is not dominated by the cen- 
tral regions where the density gradient is steep. The 
kinematic structure on large scales is also different from 
the Keplerian rotation found near the SMBH. There- 
fore, while the scenario is applicable to the narrow line 
regions of Seyferts, cannot be applicable in our case. 
The similar line width ratios of o ur results and those of 
IDe Robertis fe Osterbrockl (|1986[ ) may simply be a coin- 
cidence. 



6.3. Photo-ionization by Distributed Ionizing Sources 

The suggestion that the ionization parameter increases 
outwards naturally points to a slower decrease in the flux 
and thus to distributed ionizing sources rather than a 
central one. A number of models have been proposed 
along these lines, such as photoionization by hot evolved 
stars and by the hot X-ray emitting gas. In this section 
we first discuss the generic predictions of models with 
distributed sources, and then discuss particular models 
in more detail. 



6.2.2. Comparison to line width differences observed in 
nearby Seyferts 

Next, we consider the observed variations of line width 
between our various lines. We conclude here that the 
variations we observe are probably not related to those 
known to exist for Seyfert galaxies. 

Velocity width variations among different emission 
lines have been observed in classical nearby Seyferts and 
LINER nuclei dFilippenko fe Halpernl Il984t iFilippenkol 
[19851: IDe Robertis fc Osterbrockl 119861 : iHo et alj 11996). 
In most cases, the line widths correlate strongly with crit- 
ical density for collisonal deexcitation: lines with higher 
critical density tend to have larger widths. In a minority 
of Seyfert 2s, the line width correlates with the ionization 
potential of the ions. 

The line width differences we observe for LINERs 
is broadly consistent with that seen in local Seyferts. 
In a sample of 18 Seyfert 2 galaxies presented by 



6.3.1. Line ratio gradients from distributed sources 

Although the ionizing spectra produced by these mod- 
els differ from that produced by an AGN, the overall 
dependence of line ratios on ionization parameter and 
metallicity are very similar. Thus, as in the case of AGN, 
to explain the line ratio gradients we need the ionization 
parameter to increase outwards. Distributed sources can 
produce this trend. 

Suppose the ionizing sources are distributed like the 
stars, i.e., their luminosity density profile follows the stel- 
lar density profile. Then, we can compute the ionizing 
flux profile using the latter. Assuming that the galaxy is 
spherically symmetric, that the stellar density profile is 
p{r), and that the average number of output photoion- 
izing photons per unit time per unit stellar mass is Qo, 
the total integrated ionizing flux at distance D from the 
center of the galaxy is: 
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The integral can be solved by substituting r = D(l + e u ) 
for the r > D part and substituting r = D(l — e u ) for 
the r < D part. 

For the st ellar den s ity pr ofile, we adopt the 7-modcl 
described by iDehnenl (|1993T> : 



p(r) 



(3-7)M 



4tt 



rT(r + a) 4 "T 



(10) 



where M is the total mass, and a is a scaling factor which 
relates to the effective radius R e , depending o n the inner 
slope. 7. The 7 = 1 model correponds to the iHernquist 
( 1990) profile, the 7 = 2 model corresponds to the I Jaffe 
(1983) profile, and the 7 = 1.5 yields the best approx- 
imation of the de Vaucouleurs' i? 1 / 4 model in surface 
brightness profile. Putting this model in Eqn. [§1 we in- 
tegrate numerically to obtain the total ionizing flux as a 
function of radius for a model galaxy with an effective 
radius of 4.8 kpc and a stellar mass of 7.5 x 10 10 M Q , 
the medians for the top 25% line-emitting passive red 
galaxies in our sample. 

Figure [18] shows the resulting ionizing flux profile for 
three models with different 7 values, along with the in- 
verse square law expected from a central ionizing source. 
In the inner kpc, distributed ionizing sources will pro- 
duce a much shallower ionizing profile than the inverse 
square law. 

In Figure I19[ we divide the ionizing flux profiles by a 
gas number density profile to see how the dimensionless 
ionizing parameter will vary with radius under these dif- 
ferent models. We adopt a gas density profile that scales 

as the square root of the stell ar density profile, n g oc n~J 2 
([Mathews fc Brighcnti 2003 and references therein, also 
see discussion in £|6.2.ip , and normalize them to 100cm -3 
at 1 kpc. This warm gas density is consistent with ob- 
servations ([Heckman et al.lll989t iDonahue fc Voitlll997h 
and our assumption of pressure equilibrium between the 
warm gas (T - 10 4 K) and the hot gas (T = 10 6 - 10 7 K, 
n = 0.1 -1cm -3 ). Interestingly, for the 7 = 1.5 model, 
which gives the best fit to de Vacucouleurs' profile, the 
ionizing parameter displays the same trend as we observe, 
as shown by the [O III]/[S II] ratio profile in Fig. [H It 
not only produces the increase with radius in the central 
part but also a slow decline on the outskirts. 

6.3.2. Luminosity dependence 

A prediction of the above model is that the line 
ratio gradient s hould have a lum i nosity dependence - 
Many studies (lLauer et all 119951: |Faber et all 119971: 
Rest et all 120011: iRa vindran ath et all 120011: lLauer et al.l 
20051: iFerrarese et al.ll2006HGlass et al.ll2011[ ) have shown 
that the inner power-law slope of the stellar luminosity 
density profiles changes from —1 for bright galaxies to —2 
for faint galaxies. As shown by Figure [T9l these different 
stellar density profiles should generate different gradients 
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Fig. 18. — The integrated ionizing flux from a system of ionizing 
sources following the stellar density profile as a function of distance 
to the center of a model galaxy. The curves correspond to different 
7-models for the stellar density profile. The solid line corresponds 
to the model with the best fit to the de Vaucoulers' profile in surface 
brightness. The long dashed line corresponds to the inverse square 
law as expected in the AGN model. All models show flatter flux 
profile. 
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Fig. 19. — The ionization parameter produced by a system of 
ionizing sources following the stellar density profile in a galaxy 
shining a cloud as a function of distance from the galaxy center. 

We assumed gas density profiles of n(r) oc n 1 / 2 and normalize 
them to 100cm -3 at 1 kpc. The three curv es correspond to three 
different 7-modcls as described by Eqn. 1101 M, is the stellar mass 
of the galaxy. 
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Fig. 20. — Surface brightness profiles for the 14 early-type galax- 
ies in the Virgo cluster that would sati s ify ou r luminosity cut 
(M v < -20.4), as fit by IFerrarese et al.1 f2O0(} ). The solid gray 
curves and the dashed gray curves represent galaxies brighter and 
fainter than My of —21.3, respecti vely. The two t hick curves repre- 
sent the profiles computed for two Dolmen (1993) models, with pa- 
rameters indicated in the legend. They provide reasonable approx- 
imations to the profile difference between bright and faint early- 
type galaxies. The two vertical dotted lines indicate the range of 
apertures probed in this paper. 



in the ionization parameter. The transition point is ap- 
proximately at Mb = —20.5 (or around My = —21.3). 
Here, we investigate whether the gradients depend on 
luminosity in the expected manner. 

We divide our passive red galaxy sample at My — 
—21.3 into bright and faint samples to look for the lu- 
minosity dependence. The bright sample has a median 
stellar mass of 1.0 x 10 11 M Q and a median effective ra- 
dius of 5.8 kpc. For the faint sample, the corresponding 
numbers are 4.7 x 10 10 M(p and 3.4 kp c. 

First, we show that iDehnenl (|1993|) models with dif- 
ferent 7 values can provide reasonable approximations 
to the density profiles of bright and faint galaxies in 
our sample. To demonstrate this, we compare the 
mo dels to the sur f ace br ightness profile fits presented 
by IFerrarese et al.l (|2006l) for the 14 early-type galax- 
ies in the Virgo cluster that satisfy our luminosity cut 
(My < —20.4). In Fig. [20l the gray solid curves show 
profiles of the brighter galaxies with My < —21.3 and the 
gray dashed curves show those of the fainter ones. The 
fainter galaxies generally have steeper profiles, despite 
their s maller Sersic indices as reported bv lFerrarese et al.l 
(2006). The difference can be reasonably approximated 
by th e difference bet ween the surface brightness profiles 
of two lDehnenl (| 19931 ) models with 7 = 1.3, R e — 5.8 kpc, 
and with 7 = 1.7, R e — 3.4 kpc. 

Next, we demonstrate that the shape of integrated ion- 
ization parameter profile is not sensitive to the effective 
radius of the galaxy, but is primarily controlled by the 7 
parameter. Bright galaxies not only have shallower inner 
density profiles, but also have larger effective radii than 



Fig. 21. — Luminosity- weighted ionization parameter within 
aperture as a function of aperture radius, for five models with dif- 
ferent stellar density profiles. The density profiles are specified by 
its inner power-law slope, 7, and effective radius, R e . The gas 
density are assumed to scale as the square root of the stellar den- 
sity and is normalized to 100cm~ 3 at 1 kpc for all models. This 
figure shows that 7 controls the shape of the resulting ionization 
parameter profile while R e affects mainly the normalization. 



faint galaxies. In Figure [21] we plot models with various 
choices for 7 and R e . The ionization parameter plotted is 
the luminosity- weighted average within an aperture. For 
the luminosity weighting, we assume a power-law surface 
brightness profile with an index of —1.23, as derived from 
a fit to the [O III] profile in Fig. Combined with the 
ionization parameter profile from the model, we compute 
the luminosity-weighted average ionization parameter as 
a function of aperture radius. In all models, the gas den- 
sity is assumed to scale as the square root of the stellar 
density and they are all normalized to be 100cm -3 at 1 
kpc. 

The middle three curves in Fig.[5T]are models with the 
same R e = 4.8 kpc but different 7. They have different 
slopes in both the outer part and the inner part. They 
differ little in overall normalization. The top and bottom 
curves have the same 7 as the solid curve, but differ in 
R e . These three cases have nearly identical shapes, but 
differ significantly in their normalization. Other factors 
can impact the normalization, including Qq, the stellar 
mass, and the normalization of the gas density. Without 
knowing Qq and the normalization of the gas density, the 
normalization of the curve is free to vary. In contrast, 
only 7 controls the shape of the ionization parameter 
profile. Since we know 7 varies with luminosity, a definite 
prediction by this model is that the bright sample and 
faint sample should differ in the shape of their integrated 
[O III] /[S II] profile. 

Now we check this luminosity dependence in the data. 

4 The data show that the line emission surface brightness profile 
has only a weak dependence on the galaxy luminosity. Thus, we 
adopt the same luminosity profile for all models. 
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Fig. 22. — The aperture [O III]/[S II] ratio as a function of aper- 
ture size for the whole sample (stars), bright sample (triangles, 
M v < -21.3), and faint sample (squares, My > -21.3). The 
curves represent p redi ctions of three simple models with different 
7 parameter (Eqn. [TOt . They are not fits to the data, but a scaled 
and shifted version of the luminosity-weighted average ionization 
parameter. See text for detail. 



We divide all passive red galaxies into bright and faint 
subsamples at My = — 21.3. As in the whole sample, 
we select the top 25% galaxies in each subsample that 
have the brightest total emission line luminosity. Fig- 
ure [22] shows the [O HI]/[S II] ratios as a function of 
aperture size for the bright sample and the faint sam- 
ple separately, along with those for the whole sample. 
The curves represent the prediction of three models with 
different 7 parameter, R e and stellar mass, with the lat- 
ter two parameters adopting median values in the data. 
The models are calculated in the same way as for Fig. |2"T1 

To convert the luminosity-weighted average ionization 
parameter, we used the median stellar mass and R e for 
each model, assumed that Qq was constant, and as- 
sumed that every 0.3 dex in logf/ translates to 0.5 dex 
in [O III] / [S II] ratios. Then we shifted the models verti- 
cally by varying the gas density normalization so that the 
7 = 1.3, 1.5, and 1.7 models match roughly the normal- 
ization of the data points for the bright sample, the whole 
sample, and the faint sample, respectively. We did not 
perform an explicit fit using these data because there are 
still too many poorly-known factors in the model. Un- 
der our assumptions, the gas density at 1 kpc must be 
~ 13% lower for the faint sample than in the full sample, 
and ~ 6% higher for the bright sample. This normal- 
ization difference between the bright and faint samples 
does not have to be due to a density difference. It could 
also be due to the higher fraction of flat systems (lentic - 
ular galaxies) in the faint sample ([Bernardi et al.fe OlO). 
The stars in an intrinsically flat galaxy would be sys- 
tematically closer to the gas and yield a larger ionization 
parameter, hence higher [O III] / [S II] ratios. 

Although the normalizations match the data by de- 



FlG. 23. — Aperture [N IlJ/Ha ratios for the bright (triangles) 
and faint (squares) passive red galaxies as a function of aperture 
radius. 



sign, the shape of the models are set completely by the 
7 values. The bright galaxy sample has an increasing 
[O III] / [S II] ratio (integrated) with radius, matching the 
model prediction of a stellar density profile with a flatter 
inner slope (small 7). The faint galaxies have a much 
flatter outer slope, matching the prediction by a stellar 
profile with a steeper inner slope (large 7). In addition, 
the fainter sample displays a much steeper line ratio gra- 
dient on small radii (< 500 pc) than the brighter sample, 
matching the general trend predicted by the model. 

The data on the smallest scales do not match the data. 
This mismatch could be due to the overly simplified as- 
sumptions we made. The scaling between [O III] / [S II] 
ratio and log U may be non-linear; the inner gas density 
profile may be steeper than assumed. These could all 
change the shape of the curves. 

It is remarkable that a simple single-parameter model 
is able to predict the overall luminosity dependence of 
the [O III] / [S II] profile shape in the data. It provides 
a strong support for ionization mechanisms invoking 
sources that are distributed like the stars. 

In Fig. [221 we show the [N II] /Ha ratios as a function 
of aperture size for the bright and faint samples. The 
brighter sample always shows larger median [Nil] /Ha 
ratios than the fainter sample. This correlation be- 
tween [N II] /Ha and galaxy luminosity has been seen by 
iPhillips et al.1 (|1986| ). Our result shows that this corre- 
lation exists at all aperture scales. As we learned from 
Fig. [T71 to increase [Nil] /Ha with photoionization, we 
have to either increase the metallicity, increase the den- 
sity, or use harder ionizing spectra. The density is un- 
likely to vary by more than a factor of 10 between the 
bright galaxies and faint galaxies. And if all galaxies are 
powered by the same ionizing sources, the spectra should 
also be the same. Therefore, the difference in [N II] /Ha 
is most likely due to the gas-phase metallicity difference 
between the two samples. 
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To summarize, photoionization by distributed ioniz- 
ing sources following the stellar density profile can natu- 
rally produce the general variation of ionization param- 
eter with radius, including the sharp rise at small radius 
and gentle decline on large radius. It is also able to pro- 
duce the overall direction of the luminosity-dependence 
of the line ratio gradient. This strongly indicates that the 
spatial distribution of the true ionizing source is similar 
to the stellar distribution. 

6.3.3. Post- AG B stars 

iBinette et al.l (|1994l ) proposed that photo-ionization by 
post-AGB stars could explain the extended line emission 
in red galaxies. The spatial distribution of these post- 
AGB stars should be very similar to the overall stellar 
distribution. Therefore, based on the results of the pre- 
vious sections, the diffuse ionizing field they form can 
produce the observed line ratio gradient and its lumi- 
nosity dependence. In this section we discuss relevant 
aspects of post-AGB evolution and planetary nebulae, 
and evaluate whether they can produce enough ionizing 
photons and a sufficiently high ionization parameter. 

These stars have left the asymptotic giant branch and 
are evolving horizontally on the H-R diagram towards 
very high temperatures (~ 10 5 K) before cooling down to 
form white dwarfs. They are burning hydrogen or he- 
lium in a shell around a degenerate core. Because their 
temperatures are high enough to ionize the surrounding 
medium and plenty of material has been expelled from 
them in earlier evolutionary stages, they are often ac- 
companied by a planetary nebula. After their planetary 
nebulae disperse into the interstellar medium, the long- 
lived post-AGB stars can produce a diffuse ionizing field. 

Most of our knowledge about post-AGB stars comes 
from studies of planetary nebulae. Observations of plan- 
etary nebulae have shown their dynamical ages ar e about 
30000 years (|Schoenbernerl U9&S: iPhillipsI Tl989h . The 
time spent by stars in the post-A GB phase is a very 
strong function of their core mass (|Renzinil U983). For 
high mass post-AGB stars, the evolution is very fast: 
a post-AGB star with core mass of l. OM^ will have 
a nuclear burning time of only 25 yr (jTvlendal 119891 ) 
and fades by a factor of 10 in luminosity on a similar 
timescale. These stars evolve too fast to make their 
planetary nebulae visible for long. Thus most plane- 
tary nebulae have central stars with co re masses less 
than 0.64M Q (|Tvlenda fc Stasihskal U989). On the other 
hand, very low mass post-AGB stars (M < 0.55M Q ) 
evolve so slowly that before they raise their tempera- 
ture to 3 x 1Q 4 K (needed to ionize hydrogen) the ma- 
terial expelled in the AGB phase has completely dissi- 
pated into the in terstellar med ium. These stars, termed 
'lazy post-AGB' (|Renzinilll981h . will not appear as plan- 
etary nebulae. Therefore, considering the lifetime of 
post-AGB stars and the dynamical ages of the plan- 
etary nebulae, the central stars of planetary nebulae 
have to be post-AGB stars with core mass in a, nar- 
row range of 0.55— 0.64M^ dTvlen da fc Stasinskal ll989l 
also see lBuzzoni et al.l 120061 Fig. 15 for a nice illustra- 
tion of the mass dependent PN visibility). This core 
mass ran ge corresponds ro ughly to an initial mass of 
1— 3M Q (|Weidemaiml I2000D . Lazy post-AGB stars and 
those post-AGB stars that live longer than 30000 years 
will form a diffuse ionizing field capable of ionizing the 



neutral gas in the larger scale interstellar medium. 

However, our current understanding of the late stages 
of stellar evolution is fairly poor and we do not know 
the temperature and age distribution of these stars well. 
Most post-AGB stars observed are either hidden inside 
planetary nebulae or observed when they are not yet hot 
enough to ionize the nebula. Few hot nake d post-AGB 
stars have been obseryationally identified (jNapiwotzkil 
Il998t iBrown et al.1l200"oT: iWeston et aL|[2010h , which may 
be due to strong observational bias since they will be very 
luminous in the extreme UV but very faint in the opti- 
cal. Therefore, it is uncertain what fraction of post-AGB 
stars contribute to the large-scale photo-ionizng field and 
what fraction are hidden inside planetary nebulae. 

Fortunately, there are two clues indicating that plan- 
etary nebulae do not dominate the line luminosity in 
most of our line-emitting galaxies. First, if the line emis- 
sion is dominated by planetary nebulae, their kinematics 
should follow the stellar kinematics exactly. However, 
iSarzi et all ()2006|) showed that the ionized-gas kinemat- 
ics is decoupled from the stellar kinematics in the major- 
ity of galaxies in their sample. 

Second, we can estimate the total [O III] luminos- 
ity contributed by planetary nebulae by integrating the 
planetary nebula luminosity fun ction. We take the dou - 
ble exponential function given bv lCiardullo et all (|1989h . 

logiV(M) =0.133Af + log[l-e 3(M *~ M) ]+ const, (11) 

where M is defined as M[ ni] = — 2.5 log From] — 
13.74. The bright cut-off magnitude M* is -4.47. The 
faint cut-off magnitude is 8 mags fainter than M* 
(|Henize fc Wcstcrlund 1963). The normalization is usu- 
ally given as the total number of planetary nebulae 
within the two cut-off manitudes divided by the to- 
tal luminosity of the galaxy. We ad opt the median 
value reported by Buzzoni et al.l (|2006l ) for a sample of 
early-type galaxies, which is N = 1.65 x 10 -7 L ga i/ L & . 
With these, we found the total [O III] luminosity pro- 
duced by planetary nebulae should be L([0 III]) = 
1.35 x 10 28 L gal /L Q erg s _1 . To estimate the total PN 
light we observed through the fiber, we should use the 
fiber magnitude to derive L ga \. For the 25% passive 
red galaxies at 0.09 < z < 0.1 with the brightest to- 
tal line luminosity, the median V-band absolute magni- 
tude within the fiber aperture (derived using the fiber 
mags) is —20.09. This produces a median [O III] lumi- 
nosity of 1.21 x 10 38 erg s , much smaller than the me- 
dian [O III] luminosity among these galaxies, which is 
6.2 x 10 39 erg s _1 . Therefore, we conclude that planetary 
nebulae is a minor contributor to the total line emission 
in these galaxies. 

Can the diffuse ionizing field produced by naked post- 
AGB stars explain what we see? In this case, the sources 
are distributed like the stars and will produce the ob- 
served ionization parameter gradient and its luminosity 
dependence. The question is whether there are enough 
ionizing photons and what ionization parameter they can 
p roduce. 

IBinette et al.l (|1994h argue that the diffuse ionizing 
field has many more ionizing photons than planetary 
nebulae can produce, especially when the stellar pop- 
ulation gets older than 3 Gyr. They estimate a to- 
tal Q ~ 1 x 10 41 s^Mq 1 . The median stellar mass 
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with the fiber aperture for the 25% passive galaxies at 
0.09 < z < 0.1 with the brightest total line luminosity is 
2.1 x 1O 1O M . Assuming all post-AGB ionizing photons 
are completely absorbed and on average it takes 2.2 pho- 
toionizing photons to produce one Ha photon, this will 
produce a median Ha luminosity of 2.9 x 10 39 ergs _1 , 
which is about 1/3 of the median observed Ha luminos- 
ity 8.3 x 10 39 ergs~ 1 (before extinction correction, but 
we expect the extinction to be small). Considering the 
uncertainties on Qq and other parameters involved in 
the calculation, this can be considered as a good agree- 
ment. This luminosity is much larger than the contribu- 
tion fro m all planetary nebula e. More detailed calcula- 
tions b y IStasiriska et al.l (|2008l ) and lCid Fernandes et al.l 
(|2011| ) yield similar results. Thus, post-AGB stars can 
produce enough ionizing photons to account for most 
of the line luminosity, as long as all these photons are 
trapped inside the galaxy. This latter question is related 
with how the gas clouds are distributed relative to the 
post-AGB stars. 

In light of the decoupling of the gas kinematics from 
the stellar kinematics, let us assume that the line- 
emitting gas clouds are randomly distributed in relation 
to the post-AGB stars. In this case, we can estimate 
the ionization parameter by multiplying the Qo for post- 
AGB into Fig. [TH This yields an ionization parameter of 
\oglf = —5.2 at 1 kpc, a factor of 10 lower than what is 
re quired (log U ~ —3. 5 from Figure 1171 or -4 according 
to iBinette et al.l 119941) . Therefore, although there may 
be enough photons from post-AGB stars to produce the 
total line luminosity, in our model the light is deposited 
onto clouds far away from the ionizing sources. In con- 
sequence, the flux is significantly reduced and the result- 
ing ionization structure and line ratios are very different 
from the expectation from clouds closer to the ionizing 
sources. 

A possible solution might be that clouds closest to in- 
dividual post-AGB stars dominate in luminosity; if this 
were so, because they also have the highest ionization 
parameters, the luminosity-weighted average ionization 
parameter among all clouds would be raised. Here we 
determine that this solution is unlikely to work. To 
evaluate this possibility, we compute the average spac- 
ing between post-AGB stars using a rough number den- 
sity computed by dividing the total ionizing flux by that 
from an average star. We use an individual post-AGB 
luminosity of 10 4 L©, significantly higher than average, 
which maximizes the contribution of individual stars in 
this calculation. With this luminosity, the inter-spacing 
is around 85 pc at r = 1 kpc and increases outwards. For 
gas clouds that are randomly distributed with regard to 
the post-AGB stars, the luminosity each cloud receives 
is Flux x Area. The Flux consists of two components, 
the diffuse background Fb, and the flux from its nearest 
post-AGB star Qi/iAirr 2 ). Over a spherical volume with 
diameter equal to the inter-post- AGB spacing (r max ), the 
total luminosity due to the diffuse background is 

Air 

L bkgd = F b —r 3 max n c (A), (12) 

where n c is the number density of gas clouds and (A) is 
the average projected cloud area. The total luminosity 
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Fig. 24. — The ratio between emission line luminosity produced 
by diffuse pAGB ionizing background shining on randomly dis- 
tributed clouds and that produced by a single nearby pAGB star 
as a function of radius. This indicates that the diffuse background 
is dominant in ionizing randomly distributed clouds, confirming 
the predicted ionization parameter. 



due to an individual post-AGB star is 

Li = ( ^Anr 2 Ar)n c {A) (13) 

= Qir max n c {A) (14) 

where Q\ is the total photoionizing photon output rate 
of the star. 

Figure [24] shows the ratio L^kgd/Li as a function of 
radius for the 7 = 1.5 model. Except for the very central 
part of the galaxy, the luminosity due to the background 
is significantly larger than that due to individual nearby 
post-AGB stars. This result suggests the ionizing field 
produced by post-AGB stars is fairly smooth in most 
parts of the galaxies. The luminosity-weighted ioniza- 
tion parameter should be fairly close to what we showed 
in Fig. 1191 Granularity in the ionization field is thus 
unlikely to cause a substantially increased ionization pa- 
rameter. 

Increasing the ionization parameter requires the clouds 
to be closer to the post-AGB stars. A factor of 4 de- 
crease in average distance would probably be enough to 
bring the ionization parameter into the right ballpark. 
To achieve this, either the clouds must originate from 
progenitors of the post-AGB stars or the post-AGB stars 
must be preferentially distributed near the warm/cool 
gas. Since the gas is quite often kinematically decou- 
pled from the stars, both scenarios require that the post- 
AGBs share the same origin as the gas, rather than that 
of the main stellar population. Among all post-AGB 
stars, those with the largest core mass dominate in lu- 
minosity, which are also the youngest. If both the gas 
and the dominant post-AGB stars are associated with 
the most recent star formation episode, the dominant 
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post-AGB population might share a similar spatial dis- 
tribution and kinematics with the gas. This scenario 
would help resolve the deficit in ionization parameter. 
However, because our sample selection involves a cut in 
£>„(4000) which would exclude systems with more than 
a few percent^ of their stellar mass from a young stellar 
population (< lGyr), we consider this scenario unlikely. 
Nonetheless, it can be tested by looking at planetary neb- 
ulae kinematics to see whether they follow the stars or 
the gas. 

Another possible scenario is that the cool gas respon- 
sible for the emission indeed originates from the progeni- 
tors of the post-AGB stars. They have expanded so much 
that they no longer appear as planetary nebulae, but still 
are not as far as a randomly positioned cloud. At this 
point, they are completely dispersed in the interstellar 
medium and are carried away by the motion of the hot 
gas and thus they appear kinematically decoupled from 
the stars. In this picture, the cold gas has an internal ori- 
gin, but their kinematics is driven by the hot gas, which 
is kinematically decoupled from the stars, perhaps due 
to mergers and the collisional nature of the gas. 

Another possible solution is that the stars have a dis- 
tribution that is better resemebled by a thick disk than 
a sphere. Compared to the spherical symmetric distri- 
bution we assumed, a flatter disky distribution would 
bring the stars closer to the gas, raising the ionization pa- 
rameter. We can investigate this by looking at whether 
the stronger line emitting systems are preferentially more 
disky in morphology. We leave this for future investiga- 
tion. 

A final possible solution is that abundance of post- 
AGB stars is much larger than predicted. 

To summarize, the post-AGB star photionization 
model can naturally produce the general variation of 
ionization parameter with radius, including the sharp 
rise at small radius and gentle decline on large radius. 
It is also able to produce the overall direction of the 
luminosity-dependence of the line ratio gradient. This 
result strongly indicates that the spatial distribution of 
the true ionizing source is similar to the stellar distribu- 
tion. However, based on our current knowledge about 
post-AGB stars, the ionization parameter they produce 
would be too small, even though they may have suffi- 
cient total luminosity. The uncertainty in the number 
density of post-AGB stars is st ill too large (jBrown et al.l 
I2000ll2008t IWeston et al.ll2010f ) and observations are too 
scarce. Deeper observations and larger surveys of them 
are necessary to settle these questions. 

6.3.4. Other possible distributed photoionizing sources 

Low-mass X-ray binaries and extreme horizontal 
branch stars are two other evolved populations that 
coul d provide some ad ditional ionizing photons. How- 
ever, [SarzT^F^Ll (|2010f ) have argued that they would pro- 
duce much fewer photoionizing photons than post-AGB 
stars. Thus they are unlikely to be responsible on their 

5 The exact mass fraction of the young population tolerable by 
our D n (4000) cut depends on the assumptions used in the stellar 
population modeling. Assuming an old simple stellar population 
with solar metallicity and an age of 4.6 Gyr, which yields the ob- 
served median D n (4000) of 1.9, our D n (4000) cut at z ~ 0.1 can 
only tolerate at most 2% of its stellar mass coming from a popula- 
tion younger than 1 Gyr. 



own for the observed line emission, or to make up the 
ionization parameter deficit found above. 

Recently, high-mass X-ray binaries (HMXBs) and ul- 
traluminous X-ray sources (ULXs) (ULXs) have also 
been invoked to exp lain the LINER-like emission 
ijMcKernan et all 1201 lh . However, we do not think this 
population would solve the deficit either. First, in an 
old galaxy like those in our sample, there would be very 
few high-mass X-ray binaries, because they are only as- 
sociated with young stellar populations. Second, both 
HMXBs and ULXs are X-ray bright so they should 
have b een included in the accounting bv lEracleous et all 
(2010), who showed that extrapolating the X-ray lumi- 
nosity in the nuclear region of LINERs to the ultraviolet 
does not yield enough ionizing photons to produce the 
nuclear Ha luminosity observed. Therefore, these com- 
ponents would make at most a minor contribution. 

The hot X-ray emitting gas is also a distributed 
ionizing source that can produce LINER-like emission 
(jVoit fc Donahudll99d [Donahue fc Void 119911 ). Because 
the hot gas density approximately follows the square root 
of the stellar density, the X-ray emission should have the 
same luminosity density profile as the stars. Therefore, it 
can also produce the expected trend in line ratio gradient 
and luminosity dependence. 

However, the X-ray gas is unlikely to produce enough 
ionizing photons. For the typical galaxy in our sam- 
ple (the 25% passive red galaxy at 0.09 < z < 0.1 
with the brightest total line luminosity), with a median 
stellar mass of 7.5 x 10 10 M© and Lb = 2.8 x 10 10 L Q , 
the X-ray luminosi ty from the hot gas is on the or- 
der of 10 41 erg s _1 (jO'Sullivan et al.l 120011 ). much lower 
than the total ionizing luminosity of post-AGB stars 
(~ 10 42 erg s _1 ). Therefore, they should be subdomi- 
nant to post-AGB stars and would have an even lower 
contribution to the ionization parameter. 

6.4. Fast radiative shocks 

Shocks are prevalent in many astrophysical phenom- 
ena, such as supernova explosions, stellar winds, AGN 
jets and outflows, cloud collisions, etc. Collisional ex- 
citation in the post-shock medium can produce line ra- 
tios similar to LINERs. The fast radiative shock could 
also photoionize the precursor, unshocked regions. When 
combined with the emission lines produced in the cooling 
zone of the shock, the resulting line ra tios are similar to 
those of Seyferts ([Groves et al.ll2004bl ). In the section, 
we investigate whether the shock model can reproduce 
the trends we observed in the data. 

Because the shock-only model produces a better match 
to the LINER-like line ratios observed in these passive 
red galaxies, we only consider this model. In this model, 
the line ratios are determined by four parameters: shock 
velocity, magnetic field strength, density, and metallicity. 
Because the strong dependence of line ratios on shock ve- 
locity, if there are a wide range of shock velocities present 
in a galaxy, we should observe different widths for differ- 
ent lines in the shock-only model. The data do show 
different widths for different emission lines. The issue is 
whether the velocity dependence of line ratios can pro- 
duce the observed width difference in the right direction. 

In Section we showed that the [O III] line is on av- 
erage wider than [S II] lines by 16%. This means that 
those high-velocity line-emitting regions have a higher 
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Fig. 25. — Line ratios produced by shocks as a function of shock 
velocity for different magentic field strengths. The curves are plot- 
ted from dark black to light grey in order of increasing magnetic 
field strength. 



[O HI]/[S II] ratios than low- velocity regions. Although 
the bulk-motion velocity of clouds is not the same as the 
shock velocity, we expect in general faster moving clouds 
in a galaxy would generate faster shocks when they col- 
lide. If the lines are mostly produced in post-shock cool- 
ing zones resulting from cloud collisions, to explain the 
data would require higher [O III]/[S II] ratio to be pro- 
duced in faster shocks. 

We look at the line ratio dependence on vel ocity in the 
fast shock models given bv I Allen et al.1 (|2008| ). which are 
run using the MAPPINGS III code. Figure [25J shows 
the [O III] / [S II] ratio as a function of shock velocity for 
different magnetic field strengths. This model is run 
with solar metallicity and a pre-shock gas density of 
n = 100cm . In most part of the parameter space, 
[O III]/ [S II] decreases with increasing velocity. Only for 
B/n 1 / 2 > 10/iGcm 3 / 2 and for shock velocity between 
250 and 400km/s, does the [O III]/[S II] ratio increase 
with velocity. However, these parameter ranges do not 



produce the right velocity dependence for the other line 
ratios. The [S II] lines in our sample is narrower than 
Ha lines by 7% on average, and the [N II] lines are hav- 
ing roughly the same width as Ha lines. This requires 
the [S II] /Ha to decrease with velocity and [N II] /Ha to 
stay constant with velocity. However, the model pre- 
dicts a strongly increasing [N II] /Ha and a slightly in- 
creasing/flat [S II] /Ha with velocity in those particular 
parameter ranges, inconsistent with the requirement to 
explain the observations. Therefore, although shocks cer- 
tainly exist in most galaxies, they are probably not the 
dominant source in producing the extended line emission 
in these passive red galaxies. 

Our conclusio n agrees with the conclusion of 
iSarzi et "all (]2010( ). who also argued against the shock sce- 
nario as the dominant ionizing source based on the low 
circular velocity and low velocity dispersion observed, 
lack of morphological correlation between line emission 
structure and line ratio structure, and the flat EW dis- 
tribution. 

On the other hand, lAnnibali et al.l (|2010f ) argued that 
shocks could be important in the central regions of some 
early-type galaxies as the AGN jet-driven outflows or 
accretion onto a massive black hole could possibly reach 
the high velocities (300-500 km/s) required by the shock 
models. Shocks certainly exist in these situations, but 
having the right condition for shocks to occur does not 
mean shocks are directly responsible for the ionization 
of the gas. Further proof, such as a correlation between 
line ratio and velocity, is necessary. 

In a totally different case, in ultraluminous in- 
frared galaxies (ULIRGs), shocks could indeed be re- 
sponsible for produ cing the strong LINER-lik e emis- 
sion found there (IMonreal-Ibero et al.l 120061 120101: 
lAlonso-Herrero et al.ll201fJK As these galaxies are usu- 
ally results of major mergers, stronger and faster shocks 
are more prevalent. Star formation in these galaxies 
might also be partially responsible for the line emission. 

7. CONCLUSIONS 

In this paper, we studied the spatial distribution of 
LINER-like line emission in passive red galaxies by com- 
paring the nuclear emission luminosity measured from 
the Palomar survey with the larger aperture data from 
SDSS. We find strong evidence for line ratio gradients. 
We also find that different emission lines have differ- 
ent velocity widths, in contrast to the uniform velocity 
widths in star-forming galaxies. We have reached the 
following conclusions. 

1. In the majority of line-emitting red galaxies, the 
line emission is spatially extended and its inten- 
sity peaks at the center. The average Ha sur- 
face brightness profile can be well approximated 
by a power-law with an index of —1.28. Line- 
emitting red galaxies identified with nuclear aper- 
ture spectroscopy or those with extended aperture 
spectroscopy are essentially the same population. 

2. Line ratio gradients exist in these line-emitting 
red galaxies, with the very center having gen- 
erally larger [Nil] /Ha, [S II] /Ha, and smaller 
[O III]/[S II] than the outskirts. The [O III]/[S II] 
gradient requires an increasing ionization parame- 
ter towards larger distances. Because the cool gas 
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density is likely to fall with radius at a much slower 
rate than r -2 , an outward increasing ionization pa- 
rameter strongly disfavors AGN as the dominant 
ionizing mechanism in these galaxies. 

3. The line ratio gradient can be produced by ion- 
izing sources that are distributed like the stars. 
This model also predicts different line ratio gra- 
dient trends in bright and faint galaxies, which are 
generally matched by observations. 

4. The leading candidate for the ionizing source is 
the population of post-AGB stars. The majority 
of these stars cannot be central stars of planetary 
nebulae, but have to be naked post-AGB stars cre- 
ating a diffuse ionizing field. However, the ioniza- 
tion parameter produced by post-AGB stars falls 
short of the required value by more than a factor of 
10. Either the abundance of post-AGB stars is un- 
derpredicted or their spatial distribution has to be 
much closer to the gas clouds than assumed. The 
latter possibility would suggest a common origin of 
the gas and the post-AGB stars. 

5. Different emission lines in passive red galaxies of- 
ten have different widths. The [O III] is on average 
wider than [S II] by 16%; [N II] and Ha are wider 
than [S II] by ~ 8%. The width ratios do not vary 
as a function of aperture size. This latter result 
strongly suggests that the width ratio is not pro- 
duced by the combination of the line ratio gradient 
and rotation, but more likely due to a multiphase 
ISM in these galaxies. 

6. We considered shock models for producing these 
trends. Because line ratios produced in the cool- 
ing zone of the shocks have a strong dependence 
on the shock velocity, these models naturally pro- 
duce width differences among different lines. How- 
ever, their velocity dependence generate opposite 
width differences from what we observe. There- 
fore, shocks are strongly disfavored by our results 
as the dominant ionizing source in these passive 
red galaxies. However, it may be responsible for 
producing LINER-like emission found in ULIRGs. 

7. The systematically different [N II] /Ha ratio profiles 
between bright and faint galaxies (Fig. [23)1 sug- 
gest that the gas-phase metallicity is dependent on 
galaxy luminosity. 

Our result strongly disfavors AGN as the dominant 
ionization mechanism for the line emission in passive red 
galaxies. However, it does not mean that all LINERs 
have nothing to do with AGN. For a large fraction of 
those nuclear LINERs identified in the Palomar survey, 
accretion activity probably does exis t, as evidenced by 
the detection of compact radio core (iNagar et al J 120001) 
and X-ray point sources in their centers (|Ho et alJ l2001). 
Our result does mean that the optical line emission in 
most of them are probably powered by sources unrelated 
with AGN activity. The AGN is probably significantly 
less luminous in line emission than previously thought 
and dominates on much smaller scales than 100 pc. This 
result also helps to resolve the energy budget problem 



reported for most of these "nuclear LINERs." The large 
X-ray and radio detection fraction may be a result of 
more fuel supply in these relatively "gas-rich" early-type 
galaxies. 

For line emission found in apertures covering much 
larger area, such as in the SDSS at z > 0.02, 
the line emission is nearly always dominated by ex- 
tended emission unrelated with AGN activity. There- 
fore, most studies using line emission to derive AGN 
bolometric lumin osity for LINER -like objects using 
SPSS data fe.g. IKauffmann et all 120031: iKewlev et ali 
120061 : IKauffmann fe Heckmanl 120091: Ichoi et al.ll2009h or 
higher- z surveys (|Bongiorno et al.ll2010f) probably need 
to have their results re-inspected. The impact is proba- 
bly most significant for objects with -Lroiiii < 10 6 L Q , for 
which we have demonstrated that the ionizing sources 
are outside the nucleus. The exact threshold is also a 
function of the aperture size and galaxy luminosity. 

Although we only focused on passive red galaxies in 
our investigation, the result should also apply to LINER- 
like objects among younger red-sequence galaxies. Those 
red LINER-like objects with smaller D n (4000) probably 
also have some weak star-forming activity contributing 
to their line emission, as demonstrated in Fig. 2) 

The extended line emission is present in more than 
half of red-sequence galaxies and is much more luminous 
than most low-ionization AGNs. Based on the Palomar 
results, only the brightest few percent of low-ionization 
AGNs have a chance of detection in large aperture spec- 
troscopy data. 

Our result favors a distribution of ionizing sources that 
follows the stars, but does not confirm post-AGB stars as 
the ionizing sources. Post-AGB stars are the only source 
that has sufficient total energy to produce the observed 
emission lines. However, they fall short in the ionization 
parameter. This mystery awaits future observations to 
resolve. 

If the post-AGB stars are confirmed as the ionization 
sources, the LINERs can provide a window onto the gas 
dynamics of passive red galaxies. The total flux from 
post-AGB stars stay fairly constant with the age of the 
stellar population, except for the first Gyr after the star- 
burst. In this case, our results would then indicate that 
differing amounts of line emission in these galaxies is 
mainly an indicator of different amounts of cool/warm 
gas. We could therefore use the line strength observed 
to study the cooling and heating of warm gas in early- 
type galaxies. 

This study shows the huge amount of information we 
can learn from wide wavelength range, well-calibrated, 
high-resolution spectroscopy. It also demonstrates the 
power of large statistical samples. More detailed, spa- 
tially resolved IFU studies of nearby early-type galaxies 
are obviously the next step to confirm our results. An 
important lesson from this work is that the inclusion of 
[N II] , Ha, and [S II] in the resolved spectra was essen- 
tial to constraining the ionizing sources; indeed, in our 
case the broad wavelength range is arguably more impor- 
tant than the spatial resolution available to IFU obser- 
vations. This result motivates the use of IFU techniques 
with broad wavelength coverage to maximize the avail- 
able information. 
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